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1.  INTRODUCTION  AND  BACKGROUND 

1.1  Explosive  MHD  Generators 

Explosively  driven  MHD  generators,  as  described  In 
References  1.1  to  1.3,  operate  by  using  an  Intermediary  gas 
to  transfoxrm  explosive  energy  to  useful  electrical  output. 
Detonating  explosive  Is  used  to  shock-compress  and  energize 
argon  or  xenon  gas  to  a  dense,  high  temperature,  electrically 
conducting  plasma.  When  this  plasma  Is  passed  through  a 
magnetic  field,  electric  power  Is  generated. 

Energy  stored  In  an  explosive  Is  enormous  by  electrical 
standards.  Octol  explosive,  for  example,  co<;italns  approxi¬ 
mately  2.4  megajoules  chemical  energy  per  pound  of  explosive. 
Even  at  the  present  state  of  the  art,  and  much  remains  to 
be  done,  over  lOX  of  the  chemical  energy  Is  transferred  to 
the  Intermediary  plasma. 

1.2  Explosive  Plasma  Source  Data 

A  summary  of  recent  experimental  data  on  explosively 
driven  plasma  sources  Is  shown  In  Table  1.1.  The  data 
Includes  the  results  of  the  present  program  together  with 
results  described  In  References  1.1  to  1.3.  It  should  be 
noted  that  there  la  no  direct  method  for  measuring  the 
energy  output  of  an  explosive  driver,  and  performance  must 
be  Inferred  from  shock  wave  measurements  as  described  In 
Section  2.3  below.  Because  of  the  variety  of  plasma  pulse 
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shapes  and  the  existence  of  several  nonideal  phenomena 
occurring  In  explosive  driver  operation  (such  as  tube  expan¬ 
sion,  boundary  layers,  etc.)  the  data  In  Table  1.1  are  to 
be  considered  nominal  values  only.  Actual  values.  If  they 
could  be  measured,  might  deviate  by  as  much  as  50X  from  the 
nominal  values. 

Since  explosive  drivers  deliver  several  megajoules  of 
plasma  energy  In  a  few  tens  of  microseconds,  the  average 
theivnal  power  level  Is  very  high,  on  the  order  of  10^1  watts 
(100  gigawatts).  To  transform  this  plasma  power  to  electrical 
power  In  an  MHD  generator  requires  Intense  magnetic  fields, 
from  10  to  50  tesla.  At  these  field  levels  the  magnetic 
pressures  are  comparable  to  the  plasma  pressure.  The  Intense 
fields  are  required  only  In  the  MHD  channel,  which  Is  about 
2.5-cm  diameter  by  10-cm  long. 

1.3  Explosive  MHD  Data 

A  summary  of  experimental  data  on  the  electrical  per¬ 
formance  of  explosive  MHD  generators  is  shown  In  Table  1.2 
(References  1.1  to  1.3).  Peak  power  Is  principally  limited 
by  magnetic  field  strength  which  ranged  from  a  0.1  tesla 
permanent  magnet  to  a  3.6  tesla  self-excited  field. 

Delivered  power  varies  approximately  as  the  square  of  the 
magnetic  field. 

The  most  practical  method  for  achieving  high  values  of 
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the  magnetic  field  strength  In  single-shot  explosive  devices 
Is  self-excltatlon.  A  small  Initial  field  (.1  to  .3  tesla) 
Is  provided  from  a  permanent  magnet,  and  power  from  this 
field  Is  used  to  amplify  the  field  strength  to  full  opera¬ 
tional  values. 

1.^  Need  for  Better  Understanding  of  Nonideal  Plasmas 

Recent  self-excltatlon  experiments  (References  1.2  and 
1.3)  have  uncovered  a  serious  discrepancy  between  theory  and 
experiment  with  regard  to  plasma  resistance,  ranging  from 
a  factor  of  three  to  a  factor  of  ten.  This  discrepancy  has 
been  the  principal  technical  problem  preventing  the  achieve¬ 
ment  of  full  self-excltatlon  of  an  explosive  MHD  generator 
In  recent  experiments. 

Explosively  driven  plasmas  are  characterised  by 
extremely  high  pressures,  energies,  and  velocities  compared 
with  more  conventional  plasma  sources.  Although  the  con¬ 
ditions  are  well  suited  to  MHD  power  production,  they  are 
remarkably  difficult  for  theoretical  analysis  and  accurate 
experimental  diagnostics.  The  plasma  Is  fully  Ionized  and 
optically  dense,  so  that  most  electrical  and  optical  Instru¬ 
mentation  techniques  developed  for  thin  plasmas  are  not 
appropriate.  All  measurements  are  made  In  an  explosive  test 
environment,  with  destruction  of  most  test  components  In 
every  shot.  Theoretically,  calculation  of  explosive  plasma 
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properties  represents  a  challenge  because  of  the  density; 
free  charges  are  packed  so  close  together  that  they  resemble 
a  metallic  solid  as  much  as  a  plasma. 

1.5  Research  Program  Plan 

At  the  outset  of  this  program,  having  received  experi¬ 
mental  data  demonstrating  serious  problems  regarding  the 
prediction  of  explosively  driven  plasmas  (Reference  1.2),  It 
was  decided  to  concentrate  all  the  effort  on  achieving  a 
better  scientific  understanding  of  dense  nonideal  plasmas. 

To  this  end  a  combined  experimental  and  theoretical 
effort  was  planned.  The  experimental  program,  described  In 
Section  2  below,  concentrated  on  making  experimental  measure¬ 
ments  of  the  driver  plasma  directly  as  It  exited  the  explosive 
driver.  This  approach  eliminated  extraneous  Interactions 
with  the  MHD  channel  which  might  affect  results,  and  yielded 
the  most  well-determined  plasma  possible  from  the  explosive 
driver.  The  theoretical  program,  described  In  Section  3,  was 
concerned  with  surveying  world  literature  In  nonideal  plasma 
theory,  and  achieving  a  deeper  understanding  of  the  physical 
phenomena  relevant  to  explosive  driver  conditions. 

The  results  of  this  effort  are  encouraging,  as  discussed 
In  Section  and  It  is  probably  safe  to  say  that  quantitative 
agreement  between  experiment  and  theory  are  well  within  sight. 
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2 .  EXPERIMENTAL  PROGRAM 

2.1  Objectives  and  Experimental  Plan 

The  primary  purpose  of  the  experimental  program  was  to 
perform  diagnostic  experiments  on  high  density  plasmas  pro¬ 
duced  by  explosive  drivers.  The  principal  diagnostic  tech¬ 
niques  selected  for  this  purpose  were  an  MHD  plasma  velocity 
gage,  an  electrode-to-eleetrode  plasma  resistance  measure¬ 
ment,  an  electrodeless  conductivity  gage  utilizing  magnetic 
pickup  of  eddy  currents,  shock  wave  parameter  measurements, 
and  a  band-limited  optical  pyrometer  to  monitor  surface 
brightness  temperature  of  the  plasma. 

Three  different  types  of  explosive  drivers  were  used 
to  detennine  whether  the  mode  of  plasma  production  affected 
bulk  properties.  The  principal  concern  was  the  possibility 
of  contaminants  being  jetted  or  ablated  Into  the  plasma 
during  driver  operation.  The  three  driver  types  —  Imploding 
cylinder,  exploding  cylinder,  and  disc  —  have  such  widely 
differing  modes  of  operation  that  It  Is  unlikely  for  the 
degree  of  plasma  contamination  to  be  the  same,  and  a  compar¬ 
ison  of  measured  plasma  properties  was  expected  to  Indicate 
whether  contamination  Is  responsible  for  the  observed  low 
electrical  conductivity. 

Two  driver  gases,  argon  and  xenon,  were  selected  for 
the  test  series.  Theory  predicts  substantially  different 
temperature  and  electrical  conductivity  for  these  two  gases 
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In  an  explosive  driver,  and  this  difference  may  be  tested  by 
direct  experimental  measurement. 

The  original  experimental  plan  was  to  fire  one  shot  for 
each  of  the  two  gases  In  each  of  the  three  types  of  driver 
using  the  same  diagnostic  channel.  Constraints  of  time, 
contract  resources,  and  diagnostic  development  problems 
modified  this  plan  somewhat,  but  much  of  the  expected  data 
was  obtained.  Table  2.1  describes  the  main  features  of  the 
experiments  that  were  performed. 

2.2  Description  of  Explosive  Driver  Plasma  Sources 

Each  of  the  three  types  of  explosive  driver  —  cylindrical 
imploding,  cylindrical  exploding,  and  disc  —  transfers 
explosive  energy  to  a  dense  high  pressure  plasma  by  means  of 
a  strong  shock  wave,  but  does  so  with  differing  geometrical 
configurations  and  modes  of  plasma  processing. 

The  cylindrical  Imploding  driver  is  shown  In  Figure  2.1. 
Detonating  explosive  progressively  collapses  a  tube  containing 
the  driver  gas.  The  Imploding  tube  provides  a  cone-shaped 
gas  seal  translating  axially  at  a  velocity  equal  to  the 
detonation  velocity  of  the  explosive,  which  for  our  experi¬ 
ments  Is  8.5  km/sec  (.85  cm^sec).  The  effect  of  this  seal 
Is  to  provide  a  piston  driving  a  strong  shock  wave  into  the 
driver  gas,  which  pressurizes  and  energizes  the  plasma. 

Energy  Is  transferred  from  the  explosive  to  the  plasma  by 
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virtue  of  the  fact  that  the  Imploding  tube  must  do  work  to 
collapse  against  the  pressure  of  shocked  driver  gas,  which 
In  our  experiments  was  approximately  10  kilobars.  The 
cylindrical  Imploding  driver  represents  the  most  Ideal 
plasma  processing,  since  It  Is  essentially  a  constant  area 
shock  tube  with  constant  velocity  piston. 

The  cylindrical  exploding  driver  la  shown  In  Figure  2.2. 
It  differs  from  the  Imploding  driver  In  that  the  locations 
of  explosive  and  driver  gas  are  reversed,  and  the  detonating 
explosive  expands  a  tube  to  seal  against  an  outer  tube.  The 
moving  gas  seal  drives  a  strong  shock  Into  the  annulus  con¬ 
taining  the  driver  gas,  and  energy  Is  transferred  by  the 
expanding  tube  doing  work  against  the  shock  pressure. 
Efficiency  of  the  exploding  tube  Is  much  greater  than  the 
Imploding  tube,  but  the  gas  seal  Is  more  difficult  because 
It  Is  a  line  seal  rather  than  a  point  seal.  The  cylindrical 
exploding  driver  also  requires  a  transition  section  to 
process  the  high  pressure  plasma  from  an  annular  driver  tube 
to  a  cylindrical  MHD  channel. 

The  disc  explosive  driver  shown  In  Figure  2.3  uses  the 
exploding  driver  principle  to  drive  a  radially  converging 
shock  wave  Into  the  driver  gas.  The  resulting  plasma 
stagnates  at  the  axis  and  expands  Into  the  MHD  channel.  The 
disc  driver  Is  both  efficient  and  compact,  but  Involves  more 
complex  plasma  processing  than  either  of  the  cylindrical 
types . 
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Modes  of  potential  plasma  contamination  are  different 
for  each  of  the  driver  types.  The  Imploding  driver  has  a 
potential  jetting  problem,  similar  to  a  shaped  charge 
device,  in  which  a  Jet  of  tube  wall  material  Is  thrown 
forward  from  the  collapse  point  Into  the  plasma.  The  Jet  can 
be  minimized  and  theoretically  eliminated  by  proper  design  of 
driver  parameters,  but  there  Is  no  known  method  for  monitoring 
driver  operation  to  assure  absence  of  small  quantities  of 
Jet  contamination.  Jetting  Is  much  less  of  a  problem  in  the 
Impact  regions  of  the  cylindrical  exploding  driver  and  the 
disc  driver,  since  Jetting  is  easier  to  suppress  In  a  line 
impact  than  it  is  in  a  radially  converging  point  impact  of 
a  metal  tube. 

Ablation  from  the  driver  tube  materials  by  the  driver 
plasma  Is  another  source  of  contamination.  The  cylindrical 
Imploding  driver  is  expected  to  have  least  ablative  contami¬ 
nation  because  of  the  simple  flow  geometry  and  minimum  wall 
exposure  to  the  flowing  plasma.  The  cylindrical  exploding 
driver  Is  potentially  worse  because  of  the  transition  section, 
in  which  the  plasma  may  ablate  material  from  the  converging 
noszle  walls.  Plasma  In  the  disc  driver  first  stagnates, 
then  expands  down  an  MHD  channel,  which  may  lead  to  more 
ablation  than  would  be  expected  from  the  straight-through 
flow  of  the  cylindrical  Imploding  driver. 
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Energetic  performance  of  the  driver  and  plasma  thermo¬ 
dynamic  state  are  also  affected  by  the  driver  type.  The 
simplest  situation  Is  the  cylindrical  Imploding  driver.  In 
which  the  plasma  Is  processed  by  a  single  planar  shock  at 
constant  velocity.  The  cylindrical  exploding  driver  adds 
the  complication  of  gas  processing  and  expansion  In  the 
transition  section.  The  disc  driver  has  a  variable  velocity 
shock  In  the  disc,  and  the  plasma  Is  subsequently  stagnated 
and  reexpanded. 

2.3  Experimental  Driver  Performance 

Driver  gas  Is  energized  and  compressed  to  a  dense  plasma 
state  by  a  strong  shock  wave,  and  the  principal  means  of 
assessing  mechanical  performance  of  the  driver  and  estimating 
plasma  energy  Is  by  monitoring  the  location  and  velocity  of 
the  shock  wave. 

The  Hugonlot  shock  Jump  relationships  expressing  con¬ 
servation  of  mass,  momentum,  and  energy  across  the  shock 
front  yield  the  following  relationships  for  plasma  pressure 
and  energy; 


P=/»ouU, 


(2.1) 


(2.2) 
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Here  /i  le  Initial  loading  density  of  the  driver  gas,  u 
Is  the  plasma  flow  velocity,  and  u,  Is  the  shock  velocity. 
Plow  velocity  Is  related  to  shock  velocity  by  the  relation¬ 
ship 


u  =  (l-y)u,  (2.3) 

where  p  Is  the  shock-compressed  density  of  the  plasma. 

Shock  location  and  velocity  are  monitored  by  capped 
pins  short-circuited  by  the  shock  pressure  pulse.  Detonation 
location  Is  monitored  by  Ionization  pins  placed  against  the 
explosive,  with  ionization  phenomena  In  the  detonation  short- 
circuiting  the  pins.  The  pins  are  connected  to  a  ladder 
network  of  R-C  circuits,  so  that  short-circuiting  a  pin 
generates  a  coded  data  point  monitored  on  an  oscilloscope. 

Shota  119-17,  -18,  and  -23  were  all  Imploding  cylinder 
type  drivers  (see  Plgure  2.4  (a)).  All  the  shots  performed 
similarly,  and  an  diagram  of  Shot  119-23  Is  shown  in 
Plgure  2.5  as  an  example.  The  driver  shock  velocity  Is 
.893  cmy^sec,  about  7%  lower  than  Ideal  (.96  cm/^sec)  based 
on  measured  detonation  velocity.  It  Is  likely  that  non-ideal 
factors  such  as  tube  wall  expansion  and  boundary  layer  effects 
contribute  to  this  difference. 

Shots  119-19  and  119-20  were  both  disc  type  drivers 
(see  Plgure  2.6).  The  diagram  for  Shot  119-20  Is  shown 


-21- 


b)  Exploding  Cylinder  Type  Driver,  Shot  119-21 


Figure  2.4  Field  Photographs  of  Drivers  Used  on  Experimental  Program 
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Figure  e.5  Position- Time  Diagram  for  Imoloding  Cylinder  Type  Driver,  Shot  119-23 
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Figure  2.6  Field  Photograph  of  Disc  Type  Driver  Used  on  Shot  119-19 
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In  Figure  2.7,  indicating  the  shock  acceleration  caused  by 
radial  convergence. 

Shots  119-21  and  -22  were  both  exploding  cylinder  type 
drivers  (see  Figure  2.4(b)).  Their  x-i  diagrams,  shown  In 
Figures  2.8  and  2.9,  indicate  a  slowdown  In  growth  of  the 
shocked  plasma  region.  This  result  Is  most  likely  caused 
by  gas  seal  problems  at  the  Impact  point,  a  problem  observed 
previously  on  drivers  of  this  type.  Evidence  to  date  indi¬ 
cates  that  the  perfonnance  of  cylindrical  exploding  drivers 
Is  noticeably  dependent  on  the  circumferential  uniformity 
of  the  aluminum  expansion  tube  wall  thickness.  In  keeping 
with  this  observation  It  can  be  seen  that  the  driver  on 
Shot  119-21  performed  a  little  better  than  the  other  one. 

The  measured  Initial  wall  thickness  variation  around  the 
circumference  was  3*6J  for  this  one  as  compared  with  5-9% 
for  the  other. 

Performance  of  the  drivers  is  summarized  In  Table  2.2. 
Using  measured  shock  velocity  and  known  Initial  conditions, 
shock  pressure  and  energy  Immediately  behind  the  shock 
front  are  calculated.  These  values  typically  are  less  than 
the  nominal  value  based  on  Ideal  driver  operation,  because 
the  shock  velocity  la  typically  less  than  ideal  by  a  few 
percent,  and  both  pressure  and  energy  depend  on  the  square 
of  the  velocity. 
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Figure  2.7  Position-Time  Diagrom  for  Disc  Type  Driver,  Shot  119-20 
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Figure  2.8  Position- Time  Diagram  for  Exploding  Cylinder  Type  Driver,  Shot  119*21 
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Figure  2.9  Poeitioo-Time  Diagram  for  Exploding  Cylinder  Type  Driver,  Shot  119-22 
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2.H  Diagnostic  Channel 

A  one-lnch-dlameter  diagnostic  channel  was  used  on  all 
experiments.  Channel  material  was  Lexan  (polycarbonate 
resin  manufactured  by  General  Electric  Co.),  a  one- Inch  I.D. 
tube  with  1/8-lnch  wall  thickness.  Since  ablation  products 
affect  electrical  conductivity,  some  effort  was  Initially 
expended  In  Investigating  alternate  channel  materials. 

Glass  (silicon  dioxide)  appeared  favorable,  because  the  silicon 
ablation  products  tend  to  maintain  high  electrical  conduc¬ 
tivity,  but  practical  fabrication  problems  led  to  discontin¬ 
uation  of  the  effort.  Calculations  Indicate  that  the  Lexan 
ablation  products  reduce  plasma  conductivity  In  the  boundary 
layer  mixing  region,  and  It  was  anticipated  that  experiments 
might  tend  to  confirm  this  effect. 

On  the  shots  using  the  Imploding  cylinder  type  driver 
the  channel  was  mounted  as  an  extension  of  the  driver 
pressure  tube,  designed  to  have  a  one-inch  I.D.  The  system 
was  pressurized  to  give  a  nominal  10  kllobar  shock  pressure, 
and  consequently  the  channel  required  considerable  tamping 
to  prevent  significant  wall  motion  during  approximately 
20/1  sec  of  measurements.  To  provide  this  tamping  the 
assembled  diagnostic  channel  was  cast  In  Hydrostone  (a  high 
density  plaster,  fi  ■  1.76  g/cm3)  filled  with  09  lead  shot. 

The  mixture  was  vibrated  during  filling  to  achieve  a  high 
packing  density  of  the  shot.  This  type  of  tamper  was  used 
on  all  the  diagnostic  channels. 
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On  the  shots  using  the  disc  type  driver  the  channel  was 
mounted  to  the  driver  with  a  .001-lnch-thlck  steel  diaphragm 
Interposed.  The  channel  had  a  one-atmosphere  argon  back¬ 
ground  gas,  and  the  driver  shock  expanded  Into  this  gas.  A 
longer  duration  pressure  profile  Is  obtained  In  this  way, 
thus  slmpllflng  the  process  of  relating  measurements  made  at 
different  axial  locations. 

On  the  shots  using  the  exploding  cylinder  type  driver 
the  ring  of  shocked  gas  produced  by  the  driver  gas  was  passed 
through  a  long  (13-lnch)  tapered  section  (8  degree  taper)  and 
then  through  a  steel  diaphragm  Into  the  one-inch  diameter 
channel.  In  addition  to  using  one-atmosphere  argon  as  a 
background  gas  In  the  channel,  the  first  18  Inches  of  the 
channel  consisted  of  a  steel  barrel  to  allow  further  runout 
or  stretching  of  the  plasma  pulse. 

2.5  Design  of  Plasma  Diagnostics 
2.5.1  MHD  Velocity  Page 

The  velocity  gage  operates  as  an  open  circuit  Faraday 
MHD  generator  where  the  voltage  Is  given  by  the  equation 

V.  =Bub  (2.4) 

where  b  Is  magnetic  field  strength,  u  Is  plasma  flow  velocity, 
and  b  Is  the  Interelectrode  spacing.  With  a  known  magnetic 
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field  and  known  electrode  spacing  the  measured  voltage  provides 
the  plasma  velocity  history  by 


u(t)  = 


v,(t) 

Bb 


(2.5) 


All  diagnostic  channels  used  a  velocity  gage  with  barium 
ferrite  permanent  magnets  providing  the  magentlc  field 
(typically  B  »  0.1  tesla).  In  all  cases  a  pair  of  0.2-lnch 
diameter  brass  electrodes  protruded  0.1  Inch  Into  the  channel 
(see  Figure  2.10)  giving  a  spacing  of  0.8  Inches. 

The  electrodes  were  attached  to  a  coaxial  cable,  and 
the  output  was  viewed  across  a  50 ft  terminating  resistor 
at  the  scope.  Since  the  gas  resistance  is  ~10”2  ohm,  the 
50  ohm  termination  was  essentially  open  circuit. 

Since  the  magnets  do  not  produce  an  uniform  B-fleld 
throughout  the  volume  of  Interest  a  procedure  was  devised 
to  obtain  the  effective  B-fleld, 


(2.6) 


A  spring  driven  aluminum  bar  was  used  to  simulate  the  plasma. 
It  was  propelled  through  a  Lexan  tube  similar  to  the  channels, 
and  It  had  a  pair  of  magnets  set  up  similarly  to  the  shot 
geometry.  A  pair  of  stranded  wires  were  used  as  brush  con¬ 
tacts  for  measuring  the  voltage  across  the  bar  diameter  In 
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Figure  2.10  Brass  Velocity  Gage  Electrode  Configuration 


Figure  2.11  End  View  of  Lood  Circuit  Configuration 
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the  field.  A  lamp  and  photodiode  arrangement  was  used  to 
measure  the  bar  speed. 

The  B-fleld  on  axis  was  then  measured  using  a  Hall 
probe }  and  this  center  field  was  compared  with  the  effective 
field.  For  the  magnet  geometry  generally  used  the  ratio  was 
found  to  be 


-^=0.056.  (2.7) 

Psaatsf 

The  center  field  of  each  pair  of  magnets  on  each  shot  was 
measured  prior  to  firing,  and  the  effective  field  was 
obtained  from  this  ratio. 

2.5.2  Conductivity  via  Plasma  Resistance  Measurement 

In  addition  to  the  velocity  gage  another  Paraday  MHD 
generator  was  used  In  the  channel  to  obtain  the  plasma 
resistance  history  and  from  that  the  conductivity.  In  this 
case  the  generator  Is  provided  with  a  10  ohm  resistive 
load  which  approximately  matches  the  expected  plasma  resis¬ 
tance. 

The  voltage  appearing  across  the  load  electrodes  equals 
the  generator,  or  open  circuit,  voltage  less  the  Internal 
voltage  drop  written  as 

V,  =  Vo-iR,  (2.8) 
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where  i  Is  the  load  elroult  current  and  r.  Is  the  plasma,  or 
Internal,  resistance.  The  plasma  inductive  voltage  drop, 

L ,  Is  negligible  after  one  or  two  microseconds.  The 
plasma  resistance  history  can  be  obtained  by  measuring  , 

Vo  ,  and  )  . 

The  load  consisted  of  a  strip  of  copper  or  stainless 
steel  foil  attached  to  brass  electrodes  and  folded  to  make 
a  strip  line  to  minimize  the  load  inductance  ( 35  nH). 

Figure  2.11  shows  the  load  circuit  arrangement  Including  a 
search  coll  inductively  coupled  to  the  load  to  measure  the 
load  current.  The  load  voltage  Is  measured  across  the 
electrodes,  and  the  open  circuit  voltage  is  obtained  from 
the  velocity  gage  record.  To  do  this  account  must  be  taken 
for  the  time  lag  between  the  two  diagnostic  locations  and 
for  differences  in  B-fleld  amplitude,  the  load  circuit 
having  its  own  separate  pair  of  barium  ferrite  magnets. 

In  order  to  relate  the  plasma  resistance  to  the  net,  or 
effective,  conductivity  the  gage  factor  for  the  particular 
electrode  geometry  was  measured  in  the  lab  using  a  channel 
mock-up  with  copper  electrodes  and  a  solution  of  copper 
sulfate.  To  measure  the  conductivity  of  the  solution  the 
channel  section  was  fitted  with  copper  end  plugs.  A  current 
was  passed  through  the  entire  tube,  and  using  the  measured 
current  and  voltage  drop  as  well  as  the  length  and  cross- 
sectional  area  of  the  solution  the  conductivity  was  found  from 
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Next  a  current  was  passed  through  the  solution  via  the 
mock-up  electrodes  and  the  resulting  current  and  voltage 
drop  were  measured.  These  electrodes  were  0.4  Inches  In 
diameter  and,  like  the  velocity  gage  electrodes,  protruded 
0.1  Inch  Into  the  channel  (see  Figure  2.11).  Using  these 
measurements  and  the  measured  solution  conductivity  the 
geometry  factor,  or  gage  factor,  Is  found  as 


r  =  R<r=ya  (2.10) 

The  plasma  conductivity  on  a  shot  was  then  calculated  from 

=  •  (2.11) 


2.5.3  Conductivity  via  Eddy  Current  Measurement 

When  the  plasma  flows  through  a  gradient  In  the  B-fleld, 
there  will  likewise  be  a  gradient  In  the  MHD  voltage,  and 
this  will  result  In  circulating  currents  In  the  vicinity  of 
the  gradient  as  Illustrated  In  Figure  2.12.  A  search  coll 
placed  on  the  outside  of  the  channel  will  have  a  voltage 
Induced  In  It  by  the  variations  In  the  magnetic  field 
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produced  by  the  currents. 

Output  of  the  search  coll  Is  electronically  Integrated, 
so  that  It  directly  measures  the  field  perturbation  caused 
by  the  eddy  currents: 


V=  N AAB 

T 


(2.12) 


where  A  Is  coll  area,  n  Is  number  of  turns,  and  r  Is  the 
Integrator  time  constant.  The  magnetic  field  perturbation 
AB  Is  a  function  of  magnetic  Reynolds  number,  being  theoret¬ 
ically  linear  with  R,  at  low  values.  Thus 


AB  =  f(R.)  =  f(;4o<rub)  (2.13) 

Calibration  of  the  search  coll  In  the  lab  Is  accomplished 
using  the  spring  gun  facility  with  an  aluminum  bar  simulating 
the  plasma.  The  magpiets  and  the  search  coll  are  first  set  up 
In  the  same  geometry  as. will  be  used  In  a  shot.  T^en  the  bar 
la  propelled  through  the  mock-up  at  different  speeds  and  the 
Integrated  search  coll  output  Is  measured.  Prom  the  measured 
bar  speed  a  plot  Is  made  of  the  search  coll  output  versus  bar 
speed  (see  the  example  In  Figure  2.13). 

Since  there  are  several  orders  of  magnitude  difference 
between  the  laboratory  bar  velocities  (  ~  lo  m/sec)  and 

the  shot  plasma  flow  velocities  (  u,  10**  m/sec),  different 
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tlm«  constant  Integrators  are  used  In  the  two  cases.  Prom 
the  signal  level  obtained  on  the  shot  the  corresponding  lab 
bar  velocity  Is  obtained  from  the  plot  using, 

u.=  u,(Vj  (2.14) 


and 


Vv n  s  V,  T,  (2.15) 

where  the  l  suffixes  denote  lab  parameters  and  the  s  suffixes 
denote  shot  parameters.  The  v’s  are  Integrated  search  coll 
output  voltages,  and  the  r's  are  Integrator  time  constants. 

Since  the  same  output  signal  levels  (l.e.  Vt  )  imply 
the  same  Reynolds  number,  or  au  ,  the  plasma  conductivity 
Is  obtained  from 


The  conductivity  of  the  aluminum  bar  Is  2.49  x  10^  mho/cm, 
and  the  shot  flow  velocity  Is  obtained  from  the  velocity  gage 
or  inferred  from  the  shock  speed. 

On  Shot  119-23  a  conductivity  gage  employing  a  different 
field  geometry  was  used  In  addition  to  the  gages  described 
above.  In  this  case  a  toroidal  barium  ferrite  magnet  was 


-40- 


polarized  to  give  a  radially  Inward  B-fleld  which  In  the  plasma 
channel  results  In  an  azimuthal  MHD  voltages  and  current  flow.  The 
search  coll  was  wound  around  the  channel  and  was  located  just 
beneath  the  magnet  for  maximum  sensitivity. 

The  calibration  procedure  was  essentially  the  same  as  for 
the  other  conductivity  gage  type.  The  magnet  and  coll  were 
set  up  on  the  spring  gun  facility,  and  the  aluminum  rod  was 
used  to  simulate  the  plasma  and  thus  give  Integrated  search 
coll  output  versus  bar  speed. 

2.5.^  Optical  Pyrometer 

To  measure  the  plasma  temperature  a  small  optical  pyro¬ 
meter  was  developed  (using  Inexpensive  parts)  which  fit  onto 
the  diagnostic  channel.  The  geometrical  layout  of  the  device 
Is  shown  in  Figure  2.14.  The  device  consists  essentially  of 
an  aluminum  housing,  an  optical  train,  and  a  silicon  photo 
transistor  detector.  The  optical  train  has  a  0.0343  cm 
diameter  pinhole  drilled  in  one-mil  stainless  steel,  a  steel 
apertur,  a  90  X  width  band-pass  optical  filter  centered 
at  4500  A,  and  several  pieces  of  Pyrex  as  fillers.  In 
addition  on  one  shot  a  neutral  density  filter  was  added  as 
an  attenuator. 

The  phototransistor  type  designated  MRD  300  was  chosen 
as  the  detector  because  of  Its  low  cost,  relatively  high 
output  Into  50  ohms  ( 'v .1  volt)  and  fast  response  time 


Phototronsistor 


^ —  0.5  inch— H 


Figure  2.l4Geometrjcal  Layout  of  Optical  Pyrometer 
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(  a»  1  y«8ec)»  A  lens  on  the  front  of  the  dector  greatly 
Increases  the  effective  area  of  the  silicon  chip,  but  has  the 
disadvantage  that  the  detector  Is  extremely  angle  sensitive 
and  as  a  result  cannot  be  accurately  calibrated  In  a  straight¬ 
forward  manner.  From  a  number  of  detectors  four  were  selected 
to  be  used,  as  they  had  similar  sensitivities  and  were  the 
least  angle  sensitive  when  mounted  in  the  detector  assembly. 

Because  of  the  problems  In  achieving  reproducibility 
among  the  detector  assemblies.  It  was  decided  not  to  obtain 
absolute  calibrations  traceable  to  NBS  radiation  standards, 
but  to  proceed  with  the  manufacturers  nominal  sensitivities. 

Gage  output  may  be  thus  considered  as  a  relative  indication 
of  temperature  variation  only. 

The  manufacturer's  curves  were  used  both  for  the  trans¬ 
mission  of  the  band-pass  filter  and  for  the  spectral  sensi¬ 
tivity  of  the  MRD  300  phototransistor.  The  transmissions  of  the 
Pyrex  pieces  and  the  neutral  density  filter  at  ^500  X  were 
measured  in  the  lab.  The  pinhole  diameter  was  measured  In 
the  lab  with  an  optical  comparator. 

On  the  shots  the  phototransistor  was  provided  with  18 
volts  DC  In  series  with  a  30-ohm  output  resistor  at  the 
scope.  No  auxiliary  electronics  were  required. 
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2.5.5  Experimental  Configurations 

The  diagnostic  channel  for  Shot  119-17  (see  Figures  2.15 
and  2.16)  consisted  of  a  velocity  gage,  a  load  circuit,  a 
conductivity  gage  Just  upstream  of  the  velocity  gage,  another 
one  Just  dovmstream  of  the  load,  and  one  dummy  conductivity 
gage  further  downstream.  The  dummy  gage  consisted  of  a  search 
coll  placed  against  the  channel  with  no  magnets  to  supply  a 
field  gradient.  This  gage  was  present  to  verify  that  there 
were  no  electrostatic  or  other  stray  signals  picked  up  by  the 
search  colls.  As  the  dummy  gage  gave  no  measurable  signal. 

It  was  dropped  from  subsequent  shots. 

The  diagnostic  channel  for  Shot  119-18  was  the  same  as 
for  Shot  119-17  without  the  dummy  gage.  On  Shots  119-19  and 
-20  the  same  channel  layout  was  used  except  that  the  conduc¬ 
tivity  gage  downstream  of  the  load  was  dropped.  When  Shots 
119-17  and  -I8  were  analyzed  It  was  determined  that  the 
downstream  conductivity  gages  were  picking  up  spurious  signals 
from  the  load  circuit  currents. 

The  diagnostic  channels  for  Shots  119-21  and  -22  (see 
Figure  2.17)  consisted  of  a  velocity  gage,  two  load  circuits, 
one  with  brass  electrodes  and  one  with  tungsten  electrodes, 
a  conductivity  gage  Just  downstream  of  the  velocity  gage, 
and  an  optical  pyrometer  further  downstream.  Copper  was 
replaced  by  stainless  steel  as  the  load  material. 
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Figure  2.15  Bosic  Chonnel  Loyouf  for  Shots  119-17  thru  -20 


a)  Lexan  MHD  Channel  with  Magnets  and  Copper  Foil  Load  (Shown 
Unrolled)  Attached 


b)  MHD  Chonnel  Mounted  in  Casting  Tube  and  Ready  for  Casting 


Figure  2.16  Lab  Photogrophs  of  Typical  Diagnostic  Channel  Assembly 
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On  Shot  119-23  (see  Figure  2.l8)  the  load  circuits 
were  dropped  In  favor  of  adding  two  additional  conductivity 
gages,  one  the  usual  type  and  the  other  a  toroidal  magnet 
type.  The  magnets  for  the  velocity  gage  were  removed  to  the 
outside  of  the  casting  tube  and  the  other  square  magnets 
were  moved  out  approximately  2  cm  and  given  neoprene  buffers, 
or  shock  absorbers.  The  magnets  were  moved  thus  to  prevent 
or  delay  possible  demagnetization  by  shock  waves. 

2.6  Experimental  Results 
2.6.1  Shots  119-17  and  119-18 

The  results  of  these  two  shots  are  discussed  together 
because  of  their  similarity.  They  have  the  same  driver  type 
and  the  same  diagnostic  design  and  layout.  Experimental 
oscilloscope  traces  are  shown  In  Figures  2.19  to  2.22. 

The  expected  flow  velocity  history  at  the  velocity  gage 
location  Is  a  constant  0.85  cm//^ sec  for  approximately  15 
//sec  followed  by  a  decreasing  flow  velocity.  On  both  shots 
this  corresponds  to  16  volts  velocity  gage  output.  As  can 
be  seen  from  the  traces  there  are  departures  from  this  value 
of  up  to  4  volts  or  25f,  expeclally  In  the  first  10  ^ sec  of 
the  pulse.  Since  such  variations  In  the  gas  flow  rate  are 
unlikely,  this  Is  probably  some  spurious  signal. 

The  expected  load  voltage  and  current  for  a  constant 
flow  velocity  and  constant  plasma  resistance  are  Indicated 
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5/1  s«c 


b)  Load  Voltage  (Upper)  and  Load  Current  (Lower) 

Versus  Time 


Figure  2A9  Velocity  Gage  and  Load  Monitor  Records  for 
Shot  119-17 
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a)  Conductivity  Gage  #1  Output  Versus  Time 
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b)  Conductivity  Gage  #2  Output  Versus  Time 


Figure  £.iO  Conductivity  Gage  Records  for  Shot  119-17 
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a)  Velocity  Gage  Output  Versus  Time 
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b)  Load  Voltage  (Upper)  and  Load  Current  (Lower) 

Versus  Time 


^  * 


Figure  ?,2i  Velocity  Gage  and  Load  Monitor  Records  for 
Shot  119-18 
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In  Figure  2.23.  v.  Is  the  generator,  or  open  circuit, 
voltage,  Li  and  R,  are  the  Internal  Inductance  and  resis¬ 
tance,  and  and  are  the  external  load  Inductance  and 
resistance . 

Significant  departures  from  these  Idealized  waveforms 
occur  on  both  shots.  On  Shot  119-17  It  appears  that  due  to 
the  high  stresses  (10  kllobars  In  the  channel)  partial  shorting 
of  the  external  load  circuit  occurred  after  5 /(sec  of  flow. 

On  Shot  119-18  the  choppy  rise  In  load  volts  and  the  current 
waveform  not  returning  properly  are  Indicative  of  problems 
with  the  Integrity  of  the  load  circuit  at  these  high  stresses. 
Nevertheless  a  value  for  the  plasma  conductivity  at  5  ^sec 
was  calculated  from  these  traces  using 


a 


(2.17) 


Shot  119-17  (Argon)  gave  tr  »  40  mho/cm,  and  Shot  119-18 
(Xenon)  gave  <r  34  mho/cm.  The  latter  figure  Is  to  be 
considered  provisional . 

The  load  voltage  waveforms  may  be  subject  to  the  same 
sort  of  noise  that  appears  on  the  velocity  gage  records,  but 
In  this  case  the  variations  In  the  waveforms  may  also  reflect 
variations  In  plasma  conductivity. 

For  a  constant  flow  velocity,  plasma  conductivity,  and 
channel  size  the  conductivity  gage  outputs  should  be  constant. 
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Figure  2.23  Ideal  Load  Voltage  and  Current  Woveforms  for  Constant  Row  Velocity 
ond  Plosmo  Resistance 
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The  upstream  gages  (#1)  have  a  linear  rise  after  the  initial 
startup,  while  the  downstream  gages  have  considerably  more 
structure  including  a  sort  of  toe  or  precursor. 

The  downstream  gages  are  about  an  inch  away  from  the 
load  electrodes  which  carried  hundreds  of  amperes  on  the 
shots.  The  precursor  part  of  these  conductivity  gage  signals 
can  be  seen  to  begin  simultaneously  with  the  rise  of  current 
at  the  load.  No  such  problems  are  dlscernable  on  the  records 
of  the  upstream  gages  which  were  4  Inches  away  from  the  load 
electrodes.  Approximate  field  calculations  support  the  notion 
that  the  downstream  gage  signals  Include  a  component  due  to 
current  flowing  in  the  load  circuit. 

Since  the  conductivity  gage  outputs  are  proportional 
to  the  channel  diameter,  it  is  thought  that  channel  expansion 
is  responsible  for  the  gradual  rise  in  the  upstream  gage 
outputs.  Using  the  shock  properties  of  the  lead  and  Hydro¬ 
stone  mixture  (assuming  that  Hydrostone,  like  grout,  has  a 
5-6*  compression  at  10  kllobars)  the  channel  expansion  rate 
is  found  to  be  5.5*  in  10  >«8ec.  This  does  not  Include  con¬ 
sideration  of  the  barium  ferrite  magnets;  nevertheless  the 
Increase  in  the  conductivity  gage  signal  from  Shot  119-18 
is  about  5.2*  in  10  /<sec.  On  the  other  hand  the  Shot  119-17 
data  indicates  an  increase  of  about  29*  in  10  ji  sec .  It  is 
not  clear  how  the  channel  could  have  expanded  at  that  rate 
unless  there  were  some  voids  in  the  casting  of  lead  and 
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Hydrostone. 

Taking  the  upstream  conductivity  gage  signals  at  an  early 
time  (5/<sec  after  onset)  one  obtains  conductivities  of  50 
mho/cm  for  Shot  119-17  and  300  mho/cm  for  Shot  119-18.  The 
number  for  Shot  119-18  Is  only  an  estimate,  since  the  output 
on  the  shot  exceeded  the  range  of  the  original  calibration. 

The  estimate  was  obtained  by  setting  up  the  same  geometry  In 
the  spring  gun  with  a  different  pair  of  magnets  and  extending 
the  calibration  to  higher  bar  velocities.  Since  the  magnets 
were  slightly  different,  the  calibrations  at  a  given  bar 
velocity  are  different  (e.g.  7%  at  700  cm/sec),  and  therefore 
the  extended  calibration  Is  only  an  approximation  to  the 
original. 

2.6.2  Shots  119-19  and  119-20 

These  two  shots  both  used  disc  drivers  and  Identical 
diagnostic  channels.  The  channels  differed  from  the  previous 
ones  only  In  that  the  downstream  conductivity  gage  was 
eliminated. 

The  data  on  Shot  119-19  was  lost  due  to  early  severing 
of  the  diagnostic  cables  by  driver  debris.  A  different 
placement  of  the  cables  on  Shot  119-20  prevented  this  problem, 
and  the  scope  traces  are  shown  In  Figures  2.24  and  2.25.  The 
driver  gas  on  this  shot  was  Xenon  loaded  to  give  an  Initial 
shock  pressure  of  approximately  10  kllobars.  The  channel  gas 


b)  Load  Voltage  (Upper)  and  Load  Current  (Lower) 

Versus  Time 


Figure  2.24  Velocity  Gage  and  Load  Monitor  Records  for 
Shot  119-20 
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was  1  atmosphere  argon  separated  from  the  driver  gas  by  a 
1-mll-thlek  steel  diaphragm. 

The  expected  peak  plasma  velocity  was  about  1.1  cm//U  sec 
falling  off  to  a  third  of  that  In  about  25  ^  sec  based  on 
previous  results  for  a  similar  condlguratlon.  The  measured 
flow  velocity  history  Is  shown  In  Figure  2.26.  The  initial 
flow  velocity  was  about  1.25  cm//*  sec,  and  the  pulse  duration 
was  about  as  expected.  There  Is  again  considerable  noise  on 
the  early  part  of  the  scope  record. 

Although  the  pressure  history  In  the  channel  was  not 
measured,  It  was  a  less  severe  environment  than  In  the 
previous  shots,  since  the  driver  gas  was  expanding  Into  the 
ambient  channel  gas.  One  consequence  of  this  reduced  channel 
pressure  Is  that  the  channel  expansion  was  significantly  less, 
and  probably  no  correction  to  the  diagnostic  waveforms  Is 
necessary.  Another  consequence  was  that  the  load  circuit 
servlved  Intact  during  the  pulse.  The  load  voltage  and 
current  waveforms  show  no  evidence  of  load  Internal  breakdown. 

Figure  2,26  shows  the  plasma  conductivity  history  as 
determined  from  the  plasma  resistance.  The  plasma  resistance 
history  was  obtained  from  the  load  voltage  and  current  records. 
During  the  first  15^ sec  the  conductivity  gage  output  exceeded 
the  range  of  the  calibration.  The  two  earliest  time  points 
shown  on  the  graph  were  obtained  by  extending  the  calibration 
after  the  shot  using  a  similar  pair  of  magnets  on  the  spring 
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Conductivity  (mho/cm)  Flow  Velocity  (cin//»sac) 


a)  Measured  Flow  Velocity  Versus  Time,  Shot  119-20 


0  5  10  15  20  25  30  35 

Time  (/usec) 


b)  Measured  Plasma  Conductivity  Versus  Time,  Shot  119-20 


Figure  2.K  Plots  of  Experimental  Results  on  Shot  119-20 
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gun  facility.  The  extended  calibration  Is  credible,  since 
the  slower  bar  speed  data  matched  the  earlier  calibration 
curve.  The  point  at  6  //sec  Is  an  estimated  lower  bound, 
because  the  output  during  the  first  10 //sec  exceeded  the 
range  of  the  present  spring  gun  configuration. 

As  can  be  seen  In  the  figure  the  total  interelectrode 
resistance  at  the  load  was  decreasing  during  the  first  15 
//sec,  and  then  It  Increased  rapidly  after  that.  Further¬ 
more  the  apparent  conductivity  was  always  considerably  lower 
than  the  values  obtained  from  the  conductivity  gage. 

Some  caution  must  be  exercised  In  Interpreting  these, 
and  any  other,  results  of  the  conductivity  measurements.  If 
the  plasma  conductivity  Is  not  uniform  across  the  channel 
cross-section.  It  can  markedly  affect  the  measurements. 
Normally  there  will  be  a  boundary  layer  on  the  outer  periphery 
of  the  plasma  whose  conductivity  Is  considerably  less  than 
that  of  the  main  flow.  Unless  the  electrodes  are  so  con¬ 
figured  that  the  electrodes  tips  protrude  through  the  boundary 
layer,  the  total  resistance  between  the  electrodes  may  not  at 
all  reflect  the  conductivity  of  the  main  flow.  For  example, 
suppose  the  main  flow  conductivity  was  around  300  mho/cm 
resulting  In  a  plasma  resistance  of  approximately  2  mllll- 
ohms  for  our  electrode  geometry.  If  there  were  a  0.01  cm 
thick  boundary  layer  with  conductivity  1.  mho/cm  over  the 
surface  of  the  1  cm  diameter  electrode  (or  equivalently  a 
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0.1  cm  thick  layer  with  10  mho/cm  conductivity),  then  each 
electrode  surface  would  contribute  12.7  mllllohms  to  the 
total  resistance.  The  measured  resistance  would  then  be 
27. ^  mllllohms  which  gives  an  apparent  conductivity  of  21.5 
mho/cm . 

The  conductivity  gage  Is  also  sensitive  to  variations 
In  the  plasma  conductivity  across  the  channel  cross-section. 
Since  the  gage  search  coll  Is  placed  right  against  the 
channel  wall  for  maximum  sensitivity,  and  since  the  pickup 
Is  dependent  on  the  coupling  between  the  search  coll  and 
the  eddy  currents  In  the  plasma,  a  low  conductivity  boundary 
layer  will  result  In  a  significantly  reduced  output  signal. 

As  an  extreme  case  a  stepped  bar  was  shot  through  the 
spring  gun  to  simulate  a  sequence  of  plasmas  with  different 
zero  conductivity  boundary  layer  thicknesses.  A  scope  trace 
for  such  a  test  Is  shown  In  Figure  2.27,  together  with  a  plot 
of  output  signal  as  a  function  of  the  bar  diameter.  A  0.1 
cm  thick  boundary  layer  of  low  conductivity  will  result  In  a 
30%  reduction  In  signal,  and,  because  the  output  Is  not 
linear  with  Reynolds  number,  will  further  reduce  apparent 
conductivity.  For  example  on  Shot  119-20  the  conductivity 
gage  output  at  20  /<  sec  yielded  a  plasma  conductivity  of 
177  mho/cm  (see  Figure  2.26).  But  If  one  assumes  a  0.1  cm 
thick  boundary  layer  and  a  corresponding  30!E  signal  reduction, 
then  working  back  one  finds  that  the  main  flow  conductivity 
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a)  Conductivity  Gage  Output  Versus  Time  for  a  Stepped 
Aluminum  Bor 


Bor  Diameter  (cm) 

b)  Conductivity  Gage  Output  Versus  Bor  Diometer 


Figure  *-27  Results  of  Lab  Test  of  Conductivity  Gage  Response  Variations 
with  Bar  Diameter 
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would  actually  be  300  mho/cm. 

Both  the  plasma  resistance  measurement  and  the  eddy 
current  field  distortion  measurement  provide  spatial  averages 
of  the  conductivity  In  the  test  region.  The  details  of  this 
averaging  process  are  not  known  at  this  time,  and  await  the 
development  of  a  three-dimensional  MHD  computer  code.  In 
the  meantime  an  observation  of  the  governing  physical  mechanisms 
Indicate  that  the  plasma  resistance  gage  tends  to  weight  the 
resistive  boundary  layer  heavily,  since  resistance  of  adjacent 
layers  Is  summed.  The  conductivity  gage  measures  magnetic 
flux  from  current  flow,  and  therefore  tends  to  weight  the 
regions  of  high  conductivity.  Because  of  magnetic  field 
divergence  considerations,  the  gage  weights  current  in  the 
core  less  than  current  in  the  boundary  layer. 

2.6.3  Shots  119-21  and  119-22 

The  results  of  these  two  shots  will  be  described 
together  since  they  both  used  an  exploding  cylinder  type 
driver  and  Identical  channel  layouts.  The  only  difference 
between  them  was  the  driver  gas:  argon  for  Shot  119-21  and 
xenon  for  Shot  -19-22 .  The  drivers  were  pressurized  to  give 
a  5  kllobar  shock,  and  in  both  cases  the  channels  were  pro¬ 
vided  with  1  atmosphere  of  argon  separated  from  the  driver 
by  a  1-mil  steel  diaphragm. 
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The  scope  traces  are  shown  in  Figures  2.28  -  2.31. 

The  velocity  gage  records  were  lost,  and  after  checking  out 
the  entire  set-up  It  was  tentatively  concluded  that  the 
problem  was  an  Interralttant  failure  In  the  oscilloscope. 

For  the  purpose  of  analyzing  the  data  the  Initial  flow 
velocities  were  obtained  from  the  measured  shock  speeds  in 
the  channels  using  the  strong  shock  relation. 


y  +  t 

U.=  u 

(2.18) 

where  U,  Is  the  shock  speed,  u  is  the  flow  velocity  behind 

the  shock,  and 

y  Is  the  ratio  of  enthalpy 

to  energy .  In 

this  case  y  • 

1.23,  and  the  data  yielded 

the  following: 

Shot  # 

Us  (cm/i<sec)  u(cm//<8ec) 

Pressure 

(kilobars) 

119-21 

1.070  0.951 

1.65 

119-22 

0.974  0.865 

1.36 

The  channel  plasma  peak  pressures  are  shown  to  Illus¬ 
trate  that,  unlike  Shots  119-17  and  119-18,  the  pressure 
here  Is  low  enough  that  there  Is  negligible  channel  expan¬ 
sion  or  wall  motion.  Accordingly  the  load  voltage  and  load 
current  waveforms  show  no  indication  of  Internal  load 
shorting. 

The  #1,  or  upstream,  loads  used  brass  electrodes  of  the 
same  shape  as  on  the  previous  shots,  and  the  #2,  or  down- 
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a)  Load  ffl  Voltage  (Upper)  and  Current  (Lower) 

Versus  Time 


b)  Load  #2  Voltage  (Upper)  and  Current  (Lower) 


Versus  Time 


Figure  Z.eB  Load  Monitor  Records  for  Shot  119-21 
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a)  Load  #1  Voltage  (Upper)  and  Current  (Lower) 

Versus  Time 


SOOAmps 
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b)  Load  #2  Voltage  (Upper)  and  Current  (Lower) 


Versus  Time 


Figure  2.30  Load  Monitor  Records  for  Shot  119-22 
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a)  Conductivity  Gage  Output  Versus  Time 


J_ 

0.5 
Volts 

T 

b)  Optical  ’'yrometer  ’'hototrans later  Output  Versus  Time 


Figure  2.31  Conductivity  Gage  and  Optical  Pyrometer 
Records  for  Shot  119-22 
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stream,  loads  used  tungsten  tipped  electrodes  of  the  same 
shape.  These  extra  loads  were  Included  to  provide  compar¬ 
ative  data  on  the  dependence  of  the  Interelectrode  resis¬ 
tance  on  electrode  material  properties.  The  argument  is 
that  perhaps  ablated  electrode  material  Itself  creates  a 
high  resistance  boundary  layer  and  that  the  situation  would 
be  Improved  by  using  a  material  whose  vapor  phase  has  a  higher 
conductivity. 

The  Initial  flow  velocities  were  used  to  determine  the 
open  circuit  voltage  values  for  the  loads  and  the  flow 
velocity  values  for  the  conductivity  gages.  Lacking  complete 
velocity  histories  the  conductivities  were  determined  only 
for  5  ^  sec  after  shock  arrival.  These  results  are  tabulated 
below. 


Shot  # 

Load  01 
(Brass) 
(mho/cm) 

Load  02 
(Tungsten) 
(mho/cm) 

Conductivity 

Gage 

(mho/cm ) 

119-21 

(Argon) 

98 

68 

li<5 

119-22 

(Xenon) 

83 

51 

119 

Since  the  first  approximately  8  /x  sec  of  flow  on  both 
shots  consisted  of  compressed  argon  channel  gas  and  since  the 
driver  on  Shot  119-22  performed  measurably  poorer  of  the  two, 
the  conductivities  at  5  sec  would  be  expected  to  be  lower 
for  Shot  119-22. 


-71- 


PR-119 


Contrary  to  expectation  the  tungsten  tipped  electrodes 
did  not  show  a  higher  apparent  conductivity  than  the  brass 
electrodes.  That  the  apparent  conductivity  was  lower  may 
reflect  a  falloff  In  conductivity  with  axial  position  due 
to  progressive  contamination  since  the  tungsten  tipped 
electrodes  were  9  Inches  further  dovnstrecun  than  the  brass 
electrodes . 

On  Shot  119-22  there  Is  a  marked  Increase  In  load  voltage 
and  current  at  both  loads  and  an  Increase  In  conductivity 
gage  output  corresponding  In  each  case  to  the  arrival  of  the 
xenon  driver  plasma.  The  difference  between  the  two  shots 
can  be  Illustrated  by  comparing  conductivities  at  10  /i  sec 
obtained  from  the  conductivity  gage  outputs  and  the  Initial 
flow  velocity.  The  actual  flow  velocity  may  be  5  to  lOJJ 
lower*  so  these  results  are  only  approximate.  At  10  sec 
the  conductivity  of  the  argon  driver  plasma  (Shot  119-21)  Is 
approximately  96  mho/cm  and  the  conductivity  of  the  xenon 
driver  plasma  (Shot  119-22)  Is  approximately  l47  raho/cm. 

On  Shot  119-21  the  optical  pyrometer  phototranslster  record 
was  off  scale*  and  so  the  scope  sensitivity  was  adjusted  for 
Shot  119-22.  The  scope  record  Is  shown  In  Figure  2.31.  The 
Initial  peak  In  the  shocked  argon  channel  gas  rapidly  decays* 
and  this  Is  followed  by  an  Increased  output  upon  the  arrival 
of  the  xenon  driver  plasma.  The  rapid  fall-off  after  10/(sec 
Is  thought  to  be  due  to  shock  damage  to  the  silicon  photo- 
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transistor.  Since  the  shocked  argon  should  be  at  a  roughly 
constant  temperature,  the  decaying  signal  Is  suggestive  of 
a  growing  boundary  layer  obscuring  the  main  flow. 

One  brass  electrode  was  recovered  after  each  shot.  The 
base  of  one  tungsten-tipped  electrode  was  also  recovered, 
but  the  tip  had  been  sheared  off  and  was  not  located.  A 
brass  electrode  from  one  of  the  earlier  shots  was  also 
recovered  at  this  time.  All  of  these  electrodes  had  a 
similar  ablation  pattern  which  Is  Illustrated  In  Figure  2.32. 

Since  there  was  no  ablation  observed  on  the  remainder 
of  the  surface,  it  was  hypothesized  that  effective  thermal 
and  electrical  contact  with  the  plasma  may  have  occurred  only 
on  this  surface.  In  the  lab  sample  electrodes  were  painted 
and  a  similar  area  was  filed  away  to  expose  the  underlying 
metal.  The  gage  factor  for  these  electrodes  was  measured 
and  found  to  be  greater  than  that  of  the  fully  exposed 
electrodes  by  a  factor  of  1.7. 

If  the  hypothesis  were  true,  then  it  would  be  expected 
that  the  gage  factor  would  be  higher  yet  early  In  the  pulse 
when  only  a  small  amount  of  ablation  had  occurred.  This  Is 
consistent  with  the  observed  apparent  conductivity  for  the 
first  15  /<8ec  on  Shot  119-20  (see  Figure  2.26).  In  fact 
multiplying  the  apparent  conductivity  at  15  /wsec  by  the 
factor  1.7  gives  a  value  of  2l8  mho/cm  which  compares  well 
with  the  conductivity  gage  measurement  of  236  mho/cm. 
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2.6.4  Shot  1X9-23 

This  shot  used  an  Imploding  cylinder  type  driver  pres¬ 
surized  to  give  approximately  10  kllobar  shock  In  xenon  gas. 
The  channel  was  merely  an  Inert  continuation  of  the  driver 
pressure  tube  as  In  the  case  of  Shots  119-17  and  119-18. 

That  Is,  there  was  no  diaphragm  separating  a  channel  gas  at 
low  pressure.  The  desire  here  was  to  repeat  the  most  inter¬ 
esting  and  straightforward  configuration,  namely  xenon  plasma 
at  10  kllobars  pressure  and  0.85  cm/^u  sec  constant  flow 
velocity,  and  obtain  an  accurate  conductivity  measurement. 

The  scope  traces  are  shown  In  Figures  2.33  and  2.34. 

A  first  attempt  was  made  to  eliminate  the  noise  on  the 
velocity  gage.  The  hypothesis  was  advanced  that  stress  waves 
In  the  barium  ferrite  magnets  might  result  In  some  amount  of 
demagnetization.  Accordingly  the  magnets  were  moved  to  the 
outside  of  the  tamper  reducing  the  B-fleld  and  hence  the 
signal  level  by  a  factor  of  6  or  7.  The  magnets  for  the  con¬ 
ductivity  gages  were  moved  2  cm  away  from  the  channel  wall, 
and  pieces  of  neoprene  were  placed  In  front  of  them  to  pro¬ 
vide  additional  shock  Isolation. 

As  can  be  seen  on  the  scope  trace  (Figure  2.33)  the 
noise  level  Is  still  several  volts  while  the  expected  signal 
amplitude  Is  only  2.5  volts.  Fortunately  In  this  case  the 
flow  velocity  can  be  obtained  from  the  driver  diagnostics. 
Furthermore  It  now  appears  that  the  source  of  the  noise  Is 
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a)  Velocity  Gage  Output  Versus  Time 


b)  Optical  Pyrometer  Phototranslster  Output  Versus  Time 


Figure  S.JJ  Velocity  Gage  and  Optical  Pyrometer  Records 
for  Shot  119-23 
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a)  Conductivity  Gage  #1  (Upper)  and  Conductivity 
Gage  #2  (Lower)  Outputs  Versus  Time 


2  Volts 


b)  Conductivity  Gage  #3  Output  Versus  Time 


Figure  2,34  Conductivity  Gage  Records  for  Shot  119-23 
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something  other  than  problems  with  stress  waves  In  the 
magnets . 

Shock  waves  In  dielectric  media  commonly  give  rise  to 
transient  voltages  which.  In  the  case  of  coaxial  cables, 
result  In  noisy  signals.  Since  In  these  experiments  coaxial 
cables  extend  Into  the  high  stress  regions.  It  Is  possible 
that  such  signals  are  responsible  for  the  noise.  However, 
a  careful  mock-up  of  the  experimental  configuration  In  the 
lab  using  sharp  hammer  blows  on  confined  sections  of  coaxial 
cable  gave  only  sub-mllllvolt  noise  levels.  To  be  sure  open 
circuited  cable  gave  responses  of  several  volts,  but  with  the 
cable  end  short-circuited  to  simulate  the  channel  plasma  no 
appreciable  signals  were  observed. 

Because  the  magnets  on  the  conductivity  gages  were  moved 
away  from  the  channel,  the  B-fleld  gradient  was  significantly 
reduced.  In  order  to  get  reasonable  amplitude  signals  (  > 

0.1  volts)  It  was  necessary  to  switch  from  passive  R-C 
Integrators  to  operational  amplifier  driven  Integrators, 
particularly  for  conductivity  gage  #2  with  Its  radially 
polarized  ring  magnet. 

Although  the  Integrators  were  carefully  calibrated  and 
checked-out  In  the  lab,  apparently  some  DC  shift  at  the 
test-site  Instrumentation  bunker  resulted  In  the  outputs 
topping  out  at  2-3  volts.  As  a  result  the  first  two  con¬ 
ductivity  gages  did  not  give  a  complete  record  (see  Figure 
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2.3^(a)).  Prom  the  peak  recorded  outputs  the  lower  bounds 
on  the  conductivity  arc  3^0  mho/cm  at  conductivity  gage  HI 
and  28U  mho/cm  at  conductivity  gage  #2. 

The  downstream,  or  #3,  conductivity  gage  did  not  top 
out  probably  because  the  signal  never  exceeded  2  volts 
anyway.  The  conductivity  measured  at  5  /<  sec  from  shock 
arrival  was  130  mho/cm.  Subsequent  wall  expansion  and  flow 
velocity  changes  preclude  getting  accurate  measurements  at 
later  times. 

In  order  to  bring  the  phototranslster  output  of  the  optical 
pyrometer  down  to  the  linear  response  range  a  neutral  density 
filter  was  included  In  the  detector  optical  train  for  this 
shot.  The  response  as  shown  In  Plgure  2.33(b)  Is  more  or 
less  flat  for  approximately  3/i-<sec  and  then  rapidly  decays. 

The  onset  of  decay  corresponds  to  the  expected  arrival  of 
the  release  wave  from  the  driver  piston.  The  peak  output  is 
32  millivolts  which  corresponds  to  a  temperature  of  ^.3  eV 
based  on  the  estimated  response  curve.  The  predicted 
temperature  for  this  xenon  plasma  at  10  kllobars  Is  4.4  eV. 

2.6.5  Summary  of  Conductivity  Measurements 

The  results  of  the  conductivity  measurements  are  presented 
together  In  Table  2.3.  The  table  gives  the  results  for  each 
shot  showing  the  type  of  driver,  the  gas  used,  and  the  initial 
plasma  flow  velocity.  The  results  from  the  conductivity  gages 
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and  the  conduetlvltlea  determined  from  plasma  resistance 
measurements  are  shown  In  separate  columns.  Since  there 
appears  to  be  some  dependence  on  the  axial  location  of  the 
detector,  the  distance  of  the  detector  from  the  beginning 
of  the  Lexan  channel  la  also  shown  for  each  measurement. 

In  most  cases  conductivity  measurements  could  only  be 
accurately  made  In  the  early  part  of  the  pulse.  Accordingly 
the  conductivities  shown  In  the  table  refer  to  a  measurement  at 
5  ^ sec  after  shock  arrival  except  where  noted. 

The  most  striking  feature  of  the  data  Is  the  consis¬ 
tently  lower  apparent  conductlvltes  obtained  by  plasma 
resistance  measurements.  Some  of  the  data  supports  the 
hypothesis  that  good  electrical  contact  with  the  plasma  Is 
made  over  a  small  portion  of  the  electrode  surface  actually 
being  ablated  by  the  flow.  Then  In  terms  of  conductivity 
measurement  If  the  true  time-dependent  gage  factor  were 
known  and  applied  to  the  measured  resistance,  the  conduc¬ 
tivity  measurements  should  agree. 

Even  If  the  hypothesis  Is  true,  there  remains  the 
practical  problem  of  designing  electrodes  which  are  In  good 
electrical  contact  with  the  plasma  over  their  entire  surface. 
There  Is  not  enough  difference  In  apparent  conductivity 
between  the  brass  and  tungsten  electrodes  on  Shots  119-21 
and  119-22  to  support  the  notion  that  the  problem  Is 
material  specific.  Furthermore  previous  measurements  (Ref- 
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erenoe  1.2)  using  various  electrode  shapes  Including  wedge- 
shaped  and  finned,  or  slotted,  electrodes  likewise  gave 
apparent  conductivities  In  the  40-100  mho/cm  range.  There¬ 
fore  It  Is  not  clear  that  modifying  the  electrode  shape  will 
significantly  alter  the  performance. 

It  was  generally  observed  that  xenon  gave  higher  con¬ 
ductivities  than  argon  for  similar  conditions.  Different 
types  of  drivers  likewise  gave  significantly  different  con¬ 
ductivities.  Since  the  themodynamlc  conditions  were 
different  In  each  case,  however,  the  results  could  not  be 
correlated  as  to  driver  contamination.  In  the  theoretical 
work  It  was  found  that  for  these  plasma  conditions  the  con¬ 
ductivity  Is  a  sensitive  function  of  the  thermodynamic 
conditions.  Therefore  the  question  of  driver  plasma  contam¬ 
ination  will  have  to  be  deferred  until  each  of  the  measured 
conductivities  can  be  compared  with  theoretical  values  at 
the  appropriate  thermodynamic  conditions. 
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3.  THEORETICAL  PROGRAM 

3.1  Purpose 

The  ability  to  make  reliable  predictions  of  the  equilib¬ 
rium  thermodynamic  properties  of  explosively  driven  plasmas 
Is  needed  for  realistic  experimental  design  and  for  optimi¬ 
zation  studies.  In  the  case  of  self-excited  MHD  power  sources. 
It  Is  necessary  to  go  one  step  further  and  develop  methods  for 
determining  the  electrical  conductivity  of  the  plasma. 

The  purpose  of  the  theoretical  program  was  to  review  the 
state  of  the  art  of  predictive  methods  for  explosively  driven 
plasmas,  to  perform  preliminary  calculations  as  an  aid  to 
Identifying  plasma  parameter  regimes,  and  to  recommend  methods 
for  reducing  the  existing  discrepancy  between  experiment  and 
theory . 

3.2  Ideal  Plasma  Calculations 

The  standard  approach  In  calculating  equilibrium  plasma 
thermodynamic  state  Is  based  on  the  assumption  of  a  mixture 

of  Ideal  gases.  This  assumption  presumes  that ' particle  number 

? 

densities  are  low  enough  and  temperatures  high  enough  that 
the  average  Interparticle  Interaction  energy  Is  much  smaller 
than  the  random  (thermal)  kinetic  energy.  For  strongly 
Ionized  plasma,  Is  a  measure  of  the  average 

electrostatic  Interaction  energy  of  the  charged  particles 
and  thus  the  Ideal  gas  assumption  Implies  that  «  /  where 
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Interaction  Parameter 

y  =  e*  n'^/k  T  (3.1) 

In  Equation  3.1,  2*  la  the  mean  square  charge  and  \  Is  the 
total  charged  particle  number  density,  l.e.: 

n^*  n,  +  ^  n.  (3.2) 

“  r«€  (3.3) 

Standard  methods  for  computing  the  equilibrium  thermo¬ 
dynamics  of  mixtures  of  Ideal  gases  are  based  on  solving  a 
set  of  "Saha”  equations  together  with  charge  and  nucleil 
conservation  equations  and  are  fully  discussed  In  several 
textbooks  (e.g.  Reference  3<1)>  One  source  of  ambiguity  In 
this  procedure  occurs  In  computing  the  electronic  partition 
function,  given  by. 


9/ (3.M 

I 

where  one  has  to  devise  some  system  of  truncating  the  sum¬ 
mation  with  respect  to  energy  levels  Is  the  level 

degeneracy)  In  order  to  avoid  the  physically  unrealistic 


Infinite  sum  represented  by  Equation 

We  have  carried  out  an  Ideal-gas  calculation  for  a 
xenon  plasma  at  10  kbar  total  pressure  and  3.6  x  10^^  erg/gm 
Internal  energy  (corresponding  to  a  shock  wave  with  flow 
velocity  of  8.5  km/sec).  In  computing  these  properties  we 
assumed  that  only  the  ground  state  contributes  to  the 
electronic  partition  function  sum  In  Equation  3.4,  as  most 
of  the  truncation  schemes  drastically  curtail  the  number  of 
contributing  energy  levels  at  explosive  driver  conditions. 

Results  of  our  calculations  using  the  Ideal  gas  approxi¬ 
mation  are  shown  In  Table  3.1.  The  method  fails  a  self- 
consistency  check,  since  the  calculated  value  of  the  Inter¬ 
action  parameter  y  0.57,  well  outside  the  range  of 
applicability  of  the  theory. 

3.3  Debye-Huckel  Calculations 

The  Debye-Huckel  theory  (References  3.1  to  3,3)  has 
been  widely  used  In  computing  thermodynamic  properties  of 
weakly  nonideal  plasmas,  and  Is  the  logical  next  step  In 
our  approach.  This  method  Is  based  on  an  approximate  evalu¬ 
ation  of  the  screening  effect  of  the  surrounding  charges  on 
the  potential  of  a  given  Ion  according  to  classical  electro¬ 
statics  and  Boltzmann  statistics. 

The  screening  length  scale  or  Debye  length  of  the  plasma 


Is  defined  by 
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TABLE  3.1 

RESULTS  OF  PLASMA  CALCULATIONS  FOR  10  KBAR 
XENON  DRIVER  CONDITIONS 


Ideal  Plasma 
Approximation 

Debye-Huckel 
without  Ion 
Shielding 

Total  Pressure  (kbar) 

10.0 

10. 0 

Internal  Energy  (erg/g) 

3.60x10^^ 

3.60x1011 

Density  (g/cm^) 

.098 

.116 

Temperature  (eV) 

5.07 

M.37 

Electron  Density  (cm“3) 

7.81x1020 

1.01xl021 

Total  Ion  Density  (cm“3) 

4.<«8x1020 

5.28x1020 

Mean  Square  Charge  Z* 

1.86 

2.06 

Atom  Concentration  (cm"3) 

Xe 

2.90x10^® 

3.45x10^® 

Xe^ 

1.40x1020 

l.i<3xl020 

Xe'^+ 

2.8^1x1020 

2.94x1020 

Xe+++ 

2.31x10^9 

8.79x1019 

Xe+++^ 

1.06x1018 

2.94x1018 

Interaction  Parameter  y  *565  .783 
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=  [  kT/4nc*  (n^  ^  f  Bt  n,)j  (3.5) 

Electrostatic  Interactions  led  to  a  lowering  of  the  Ionization 
potential  given,  for  an  Ion  of  charge  ,  by 

A  I,.  =(z^+i)  c*/Dc  (3.6) 

and  a  reduction  of  the  total  pressure 

A  P  ^  c’  C  n,  z]  Hi  ) / 6  Dc  (3.7) 

Thus  In  rolvlng  for  the  thermodynamic  state  using  the 
Debye-Huckel  procedure,  one  essentially  repeats  the  steps 
taken  In  the  Ideal  gas  approximation  approach  but  with  new 
Ionization  potentials  1^'**  -  I,.  -  and  p"'"  =  -  a  p  . 

In  calculating  the  electronic  partition  functions,  one 
computes  the  sum  In  Equation  3.4  to  ^  ,  thus  removing 

the  ambiguities  associated  with  the  truncation  of  the  Infinite 
sum. 

In  attempting  to  calculate  the  thermodynamic  properties 
of  a  xenon  plasma  with  10  kbar  total  pressure  and  3.6l  x  10^^ 
ergs/gm  Internal  energy,  we  encountered  thermodynamic  Insta¬ 
bilities  (characterized  by  calculated  total  pressures 
becoming  negative).  Since  the  question  of  full  participation 
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by  the  ions  In  the  dynamic  coulomb  shleldlnp;  of  a  plasma  Is 
subject  to  some  controversy  (Reference  3.34),  we  decided  to 
try  the  Debye-Huckel  calculations  assuming  screening  by 
electrons  only  le,  with  the  Debye  length  taken  as 

Dg  =  [  kT/4- rr  e.*  HgJ  (3.8) 

This  approach  leads  to  thermodynamically  stable  results  which 
are  also  shown  in  Table  3.1. 

The  question  of  ion  participation  in  the  electrostatic 
shielding  was  considered  by  Devoto  (Reference  3.34)  who 
computed  electrical  and  thermal  conductivities  of  argon  at 
1  atm  pressure  by  solving  the  Boltzmann  collision  term  with 
the  shielded  coulomb  potential  using  both  definitions  of  the 
Debye  length.  Comparisons  with  calculations  of  Williams  and 
Dewitt  (Reference  3.28),  who  used  a  quantum  mechanical 
formulation  in  which  the  dynamic  shielding  process  was 
treated  correctly,  revealed  better  agreement  if  shielding  by 
ions  was  ignored. 

The  Debye-Huckel  technique  is  generally  expected  to  be 
valid  for  high  temperatures  and  low  densities  (l.e.  a  weakly 
nonideal  gas)  or,  equivalently,  in  the  limit  that  the  Debye 
sphere  contains  a  large  number  of  charged  particles.  A  less 
restrictive  condition  for  its  applicability  has  been  derived 
(Reference  3.2)  from  a  classical  statistical  mechanical 
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study  of  weak  electrolyte  solutions.  This  condition  can  be 
expressed  In  the  form, 


He  •  >^e  ^  ^  Cs  ^  Dc  ) 


(3.9a) 


or,  equivalently, 

s  ot  (n,  Zn-)  ^  j  (3.9b) 

Y  ^  =0.47  (3.9c) 

This  restriction  Implies  that  on  the  average  it  Is 
necessary  to  have  at  least  one-sixth  of  a  charged  particle  In 
the  Debye  sphere  for  the  Debye-Huckel  theory  to  be  valid.  It 
needs  to  be  pointed  out,  however,  that  both  the  Debye-Huckel 
theory  and  the  condition  3.9a,  b,  c  for  Its  applicability 
are  based  on  the  assumption  of  a  classical  mechanical  descrip¬ 
tion  of  charged  particle  Interaction.  As  discussed  In  later 
sections,  quantum  mechanical  effects  play  an  important  role 
In  a  strongly  nonideal  gas  and  clearly  the  "real"  region  of 
applicability  of  the  Debye-Huckel  theory  Is  reduced  to  values 
of  Y  where  quantum  mechanical  effects  begin  to  become 
appreciable. 
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Using  the  Debye-Huckel  theory  with  electronic  partici¬ 
pation  in  the  shielding  process,  the  thermodynamic  parameters 
of  Table  3.1  were  calculated.  The  number  of  electrons  in  a 
Debye  sphere  is  about  0.5,  which  is  within  the  range  of 
applicability  of  the  theory.  The  interaction  parameter  y 
is  equal  to  0.79,  indicating  a  strongly  nonideal  plasma. 

The  Debye-Huckel  expression  has  been  shown  to  be  the 
first  order  term  in  a  perturbation  solution  of  the  classical 
statistical  mechanical  formulation  for  nonideal  plasmas 
(Reference  3*2)  and  much  work  has  been  done  to  evaluate 
higher  order  terms  (see  for  example.  Reference  3.39).  How¬ 
ever,  this  approach  is  not  useful  for  the  strongly  nonideal 
plasmas  of  Interest  to  us  since  any  Improvement  on  the  Debye- 
Huckel  approach  must  take  account  of  the  quantum  mechanical 
effects . 

The  quantum  mechanical  effects  can  be  generally  sub¬ 
divided  into  two  types;  electron  degeneracy  effects,  and 
effects  associated  with  the  quantum-mechanical  phenomena 
associated  with  strong  (l.e.  short-range)  charged  particle 
collisions.  One  approach,  widely  used  in  reentry  physics 
and  nuclear  fireball  calculations  about  a  decade  ago,  is 
based  on  the  ion-sphere  theory  of  plasma  interaction.  In 
this  approach,  one  accounts  for  electron  quantum  degeneracy 
effects  by  treating  the  electrons  within  the  framework  of  a 
finite  temperature  Thomas-Perml  statistical  model.  In 
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fact  a  general  theory  (Reference  3 >6)  Is  available  which 
provides  the  Debye-Huckel  and  the  Ion-sphere  solutions  as 
limiting  cases  (high  temperature,  low  density  limit  for  the 
former  and  low  temperature,  high  density  limit  for  the  latter). 
However,  this  approach  Is  Inadequate  for  plasma  conditions  of 
Interest  to  us  where  It  Is  essential  to  account  for  quantum 
mechanical  collision  effects.  Indeed  the  breakdown  of 
thermodynamic  stability  associated  with  a  Debye-Huckel  theory 
based  on  full  Ion  participation  In  the  screening  process  is 
an  Indication  of  the  need  to  account  for  the  quantum  colll- 
slonal  effects. 

3.4  Parameters  Characterizing  Plasma  Regime 

It  Is  useful  at  this  stage  to  consider  the  parameters 

that  characterize  the  plasma  regime.  We  consider  a  multiply 

ionized  plasma  for  which  +  Z  is  the  total 

< 

charged  particle  number  density  and  z*  =  Z  2  -  n^)/n^ 

is  the  mean  square  charge  on  a  charged  particle.  There  are 
five  primary  length  scales  as  follows: 

Debye  length: 

(kT/-4  TT  e*  n^)'^  =  7.34-3x10'“  (T/z*n^)'^  {'/i )  (3-10) 
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Average  Impact  Parameter: 

6.  *  i*  cV3kT  =  4.8  00  P/t  (a) 

Thermal  DeBroglle  Wavelength  of  Electrons; 

h /(  kT)'^^  ^  L952/t'^(A) 

Average  Charged  Particle  Separation 

Average  Electron  Separation 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


The  Debye  length  Is  a  measure  of  the  effective  range  of 
the  coulomb  potential.  As  mentioned  In  the  previous  section, 

It  can  be  defined  In  terms  of  electron  participation  In  the 
shielding  process  (  )  or  participation  by  both  electrons 

and  Ions  (  D-D^  ).  fa.  Is  the  average  classical  Impact 
parameter  for  90®  scattering  and  qualitatively  delineates 
strong  and  weak  scattering  regions. 

Prom  these  primary  length  scales  we  construct  the  following 
Important  dimensionless  parameters  with  which  we  can  characterize 
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the  plasma  regime. 

Interaction  Parameter: 

Y  =  z*  e*  -  1.440x10-'  (3.15) 


Plasma  Parameter 

a  a:  =  ^  y-^  (3.16) 

Number  of  particles  in  Debye  Sphere 

^  Di  =  A,  /g  =  TO  r'*'^  (3.17) 

Quantum  Diffraction  Parameter 

^e/D^  (3-18) 

Electron  Degeneracy  Parameter 

^eAe  (3.19) 

A  classical  nonideal  plasma  is  completely  characterized 
by  the  interaction  parameter  y  •  Both  the  plasma  parameter 
A  and  the  number  of  charged  particles  in  the  Debye  sphere, 

Nq  can  be  expressed  in  terms  of  it.  The  importance  of 
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quantum  mechanical  effects  is  characterized  by  two  additional 
parameters,  representing  electron  degeneracy  effects, 

and  the  quantum  diffraction  parameter  X^/d  which  charac¬ 
terizes  the  quantum  nature  of  the  colllslonal  Interaction  of 
electrons  with  other  electrons  or  ions.  It  is  useful  to  note 
that,  in  the  weakly  nonideal  plasma  limit,  the  five  primary 
length  scales  are  ordered  according  to: 

(3-20) 

and  hence  in  this  regime  (  |^«  1  ,  X./Pc  .  «  I  )  a 

theoretical  formulation  based  on  classical  mechanical  descrip¬ 
tion  of  charged-particle  scattering  and  neglect  of  degeneracy 
is  adequate. 

The  xenon  plasma  calculations  are  a  useful  guide  in 
characterizing  the  plasma  regimes  encountered  in  typical 
explosive-driven  MHD  power  sources.  The  length  scales  and 
dimensionless  parameters  have  been  evaluated  on  the  basis 
of  calculated  equilibrium  thermodynamic  properties  of  a  10 
kbar,  8.5  km/xec  xenon  plasma,  and  are  shown  in  Tables  3.2 
and  3.3*  We  have  computed  the  Debye  length  (D  ),  plasma 
parameter  (A)  and  number  of  charged  particles  in  the  Debye 
sphere  {  )  in  two  ways:  accounting  for  both  ions  and 

electrons  (subscript  "c")  and  accounting  for  the  electrons 
only  (subscript  "e"),  and  this  would  identify  the  differences 
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TABLE  3.2 

CALCULATED  PLASMA  LENGTH  SCALES  POH  10  KBAR 


XENON  DRIVER 

CONDITIONS 

Ideal  Plasma 

Debye  Huckel 

Debye  Length  with 

Ion  Shielding,  A 

3.50 

2.76 

Debye  Length  without 

Ion  Shielding,  % 

5.99 

4.90 

Electron  DeBroelle 
Wavelength,  A 

.867 

.933 

Electron  Separation,  X 

10.9 

9.98 

O 

Impact  Parameter,  A 

1.76 

2.26 
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TABLE  3.3 

CALCULATED  DIMENSIONLESS  PLASMA  PARAMETERS 
FOR  10  KBAR  XENON  DRIVER  CONDITIONS 


Ideal  Plasma 

Debye-Huckel 

Interaction  Parameter  y 

.565 

.783 

Plasma  Parameter 

With  Ion  Shielding 

1.99 

1.22 

Without  Ion  Shielding 

3.^1 

2.17 

Charged  Particles  in 

Debye  Sphere 

Ions  and  Electrons  A/- 

.222 

.136 

Electrons  Only  A/_ 

.168 

.496 

Degeneracy  Parameter 

.0798 

.0935 

Quantum  Scattering  Parameter 

.2A8 

.338 

\./0, 
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resulting  from  the  choice  of  shielding  model  used. 

The  results  indicate  that  the  plasma  is  strongly  nonideal, 
with  an  Important  a.uantum  mechanical  contribution  to  charged- 
particle  scattering  and  with  a  significant  electron  degeneracy 
effect.  It  is  quite  clear  that  the  use  of  any  theory  such  as 
the  Debye-Huckel  or  the  Ion-sphere  (or  an  interpolation 
between  the  two)  that  does  not  account  for  both  types  of 
quantum  mechanical  effects  will  not  yield  a  satisfactory 
description  of  the  equilibrium  thermodynamic  properties  of 
the  plasmas  occurlng  in  explosively-driven  MHD  power  sources. 

3.5  Survey  of  Current  Research  in  Nonideal  Plasmas 

In  this  section,  we  will  summarize  developments  in 
research  on  equilibrium  thermodynamic  properties  and  elec¬ 
trical  conductivity  of  nonideal  plasmas  published  during  the 
last  decade.  A  selected  number  of  prior  efforts  will  be 
also  discussed  in  order  to  place  the  more  recent  work  in 
perspective. 

3.5.1  Equilibrium  Thermodynamic  Properties 

We  have  tabulated  a  list  of  publications  relating  to 
the  determination  of  equilibrium  thermodynamic  properties 
of  nonideal  plasmas  in  Table  3.^>  corresponding  to  References 
3.2  to  3.23.  Perhaps  the  most  Important  contribution  among 
these  is  the  exhaustive  review  of  strongly  nonideal  plasmas 
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by  Norman  and  Starostin  (Reference  3.1*»).  This  review, 
published  In  1970,  covered  state  of  the  art  thermodynamic 
property  calculation  techniques  for  both  weakly  nonideal 
(  /■«l  )  and  strongly  nonideal  plasmas  (  y'  ~  Various 

questions  in  nonideal  plasma  theory  are  discussed,  Including 
the  Inadequacy  of  classical  models,  quantum  diffraction 
effects,  electron  degeneracy,  phase  transitions  and  thermo¬ 
dynamic  stability. 

It  is  noted  that  classical  approaches  such  as  the  Debye- 
Huckel  based  on  classical  statistical  mechanics  (Reference 
3.1),  become  increasingly  Inaccurate  as  the  Interaction 
parameter, y-  ,  increases  beyond  about  0.1.  In  this  regime 
(which  Includes  explosively  driven  plasmas)  It  Is  necessary 
to  take  Into  account  quantum  mechanical  effects,  such  as 
diffraction  and  electron  exchange  effects  In  charged  particle 
collisions  and,  when  the  electron  number  density  Is  suffi¬ 
ciently  high,  electron  degeneracy  effects.  Only  by  appro¬ 
priately  accounting  for  such  quantum  mechanical  effects  can 
thermodynamic  stability  be  guaranteed  and  phenomena  such  as 
phase  transitions  and  plasma  metallization  be  correctly 
treated . 

It  Is  noteworthy  that  prior  to  the  late  1960's  the 
Important  advances  In  research  on  thermodynamics  of  nonideal 
plasmas  was  predominantly  of  U.  S.  origin.  The  emphasis 
during  this  period  was  on  weakly  nonideal  plasmas  relevant 
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to  re-entry  type  of  plasma  environment.  In  more  recent 
years  the  situation  has  changed  remarkably  and  In  contrast 
to  the  paucity  of  effort  In  the  U.  S.  and  other  western 
nations,  the  Soviet  Union  appears  to  be  engaged  In  a  deter¬ 
mined  effort  to  understand  strongly  nonideal  plasmas. 

This  vigorous  Soviet  effort  has  been  directed  towards 
both  theoretical  and  experimental  studies  and  a  concomitant 
emphasis  on  the  derivation  of  useful  seml-emplrlcal  results. 
The  Soviet  theoretical  work  Is  being  pursued  primarily  by 
Norman's  group,  and  It  Is  useful  to  summarize  the  status  of 
their  theoretical  effort: 

1.  Because  of  the  lack  of  a  suitable  small  parameter 
In  the  nonideal  plasma  regime.  It  is  necessary  to  rely  on 
involved  numerical  approaches  such  as  the  Monte  Carlo  method 
for  a  theoretical  solution  of  the  quantum  statistical  mechan¬ 
ical  description  of  strongly  nonideal  plasmas. 

2.  The  fundamental  problem  Is  the  choice  of  charged- 
particle  Interaction  laws  that  effectively  account  for  quantum 
effects  In  a  nonideal  plasma.  The  Soviets  have  devoted  a 
great  deal  of  effort  to  devising  a  "pseudopotential"  which 
correctly  accounts  for  quantum  mechanical  effects  in  short- 
range  collisions  but  at  the  same  time  reproduces  the  coulomb 
potential  at  large  distances. 
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3.  Most  of  the  theoretical  studies  so  far  have  been 
for  model  plasmas  without  taking  account  of  degeneracy.  The 
plasma  model  most  commonly  chosen  Is  one  In  which  the  atoms 
are  treated  separately  as  an  Ideal  gas  co-exlstlng  In  lonl- 
zatlonal  equilibrium  with  a  nonideal  subsystem  composed  of 
the  Ionic  constituents  and  free  (continuous  spectrum)  electrons. 
The  Interaction  of  the  Ions  with  the  continuous-spectrum 
electrons  Is  described  by  means  of  the  pseudopotentlal . 

*♦.  Using  the  Monte  Carlo  method,  the  equation  of  state 
and  Ionization  equilibrium  of  the  model  plasma  are  computed 
over  a  wide  range  of  x  •  Results  from  experimental  studies 
(e.g.  Reference  3»18)  of  nonideal  plasmas  are  being  used  to 
estimate  undetermined  parameters  In  the  pseudopotential  and 
to  aid  in  derivation  of  similarity  laws. 

One  Important  U.  S.  contribution  In  recent  years  Is  that 
by  Rogers  and  DeWltt  (Reference  3*20)  who  have  developed  a 
quantum  statistical-mechanical  formulation  for  the  equations 
of  state  of  coulomb  gases.  However,  this  approach  requires 
extensive  numerical  computations  and,  within  the  framework 
of  a  number  of  severe  approximations,  has  so  far  only  been 
applied  to  the  determination  of  the  equation  of  state  of 
gaseous  hydrogen.  It  Is  interesting  to  note  that  in  the 
hydrogen  plasma  calculations,  the  Debye-Huckel  relations 
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are  found  to  be  adequate  for  densities  not  exceeding 
10"^  gm/cm^. 

In  conjunction  with  a  strong  theoretical  effort,  the 
Soviets  have  and  are  continuing  a  vigorous  experimental 
program.  As  mentioned  above,  experimental  data  is  needed 
for  determination  of  undetermined  parameters  in  the  theoret¬ 
ical  model.  Experimental  data  is  also  valuable  in  developing 
seml-emplrlcal  expressions. 

A  great  deal  of  the  experimental  effort  has  concerned 
cesium  because  of  the  relative  ease  with  which  strongly  non- 
ideal  plasma  conditions  can  be  obtained.  Fortov  (Reference 
3.17)  has  considered  a  large  number  of  existing  methods  for 
computing  the  thermodynamic  properties  of  a  nonideal  cesium 
plasma  and  computed  these  properties  for  pressures  and  tem¬ 
peratures  In  the  range  5-150  atm.  and  3,000  to  15,000®  K 
(the  largest  value  of  ^  Is  1.08  In  this  range!  His  results 
Indicate  striking  differences  In  the  results  obtained  by  the 
different  theories. 

This  result  Illustrates  the  Importance  of  satisfying 
oneself  that  the  approach  used  In  thermodynamic  calculations 
is  appropriate  to  the  plasma  regime  of  Interest,  particularly 
If  one  Is  Interested  In  reliable  computations  for  the  elec¬ 
trical  conductivity  which  Is  strongly  dependent  on  thermo¬ 
dynamic  properties  such  as  temperature,  electron  density  and 
average  Ionic  charge. 
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Ermokhln  (Reference  3*15)  has  conducted  experiments  with 
cesium  under  conditions  such  that  ,  the  number  of  particles 
In  the  Pebye  sphere  varies  from  about  0.6  to  ^<50.  He  has 
shown  that  the  effect  of  coulomblc  Interaction  In  this  regime 
Is  such  that  a  fairly  reliable  description  Is  obtainable  If 
one  Interpolates  between  the  weakly  nonideal  Chapman-Enskog 
approach  for  large  and  the  theory  of  electron  Interactions 
with  Impurity  Ions  In  a  semiconductor  (Reference  3-^0)  valid 
for  «  /  . 

Experiments  with  Inert  gases  such  as  xenon  and  argon 
which  are  of  direct  Interest  to  us,  have  been  comparatively 
limited  and  usually  not  under  the  conditions  of  very  strong 
nonidealities  explored  In  cesium  plasmas.  This  Is  probably 
because  of  the  experimental  difficulties  In  producing  strong 
nonidealities  In  these  gases. 

3.5.2  Electrical  Conductivity 

A  review  of  the  literature  on  methods  for  calculating 
electrical  conductivity  of  a  fully  Ionized  nonideal  plasma 
Is  given  In  Table  3.5»  corresponding  to  References  3.2^ 
to  3.  38. 

Historically,  rigorous  kinetic  theory  approaches  to  the 
calculation  of  the  electrical  conductivity  have  been  based 
on  the  Chapman-Enskog  type  of  approximations  to  the  collision 
Integral  In  the  Boltzman  equation.  This  approach  works  well 
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SUMMARY  OF  RESEARCH  ON  ELECTRICAL  CONDUCTIVITY  OF  NONIDEAL  PLASMAS 
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1975  Mltin  USSR  Measurement  of  conductivity  of  argon  and  copper  In  a 

pulsed  arc  at  conditions  corresponding  to  1.6 
Finds  Spltzer  theory  In  agreement  with  measurements. 
Develops  ad  hoc  seml-emplrlcal  improvement  of  Spltzer 
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for  dilute  neutral  pases  and  can  be  applied  to  plasmas 
provldlnF  some  sort  of  shielding  Is  applied  to  the  coulomb 
potential  to  prevent  the  long  range  divergence  Induced  by  the 
pure  coulomb  potential.  Usually  the  potential  Is  "cut  off" 
at  the  Debye  length  or  one  uses  a  shielded  coulomb  potential 
In  which  the  Debye  length  Is  the  characteristic  shielding 
distance.  An  equivalent  approach  is  to  consider  the  Fokker- 
Planck  equation.  For  fully  Ionized  plasmas,  the  classic 
electrical  conductivity  expression  Is  that  of  Spltzer-Harm  (Ref 
erence  3.2^)  which  provides  a  simple  expression  valid  pro¬ 
vided  A  A  »  I  .  There  have  been  several  efforts  to  Increase 
the  range  of  applicability  of  this  approach  based  on  both  the 
classical  Boltzman  equation  technique  (Reference  3*^1)  and  on 
quantum  mechanical  plasma  transport  theory  (Reference  3.28). 
However,  none  of  these  approaches  extend  the  region  of  appli¬ 
cability  to  better  than  A  s  3  and  are  thus  not  appropriate 
for  strongly  nonideal  plasmas. 

The  problem  arises  from  the  Inadequacy  of  the  classical 
collision  models  used  by  all  these  approaches  since,  under  con¬ 
ditions  of  stronge  nonideality,  quantum  mechanical  short  range 
collision  effects  make  an  Important  contribution  to  the 
transport  integrals.  An  approximate  quantum  mechanical 
approach  has  been  developed  by  Rogoff  (Reference  3.30)  who 
uses  a  Born  approximation  with  the  Debye-shielded  potential 
within  the  framework  of  a  fourth-order  Chapman-Eskog  approx- 
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Imatlon  for  a  fully  ionized  plasma.  Moreover,  his  results 
are  expressible  In  terms  of  a  simple  analytical  expression 
In  the  high  temperature  limit.  This  result  is  clearly  of 
great  value  under  conditions  of  small  values  of  A  where  the 
classical  collision  techniques  based  on  the  Boltzman  equation 
break  down  dramatically.  However  we  do  not  yet  have  a  reliable 
indication  of  the  accuracy  of  this  method. 

A  more  accurate  theory  of  electrical  conductivity  for 
nonideal  plasmas  has  been  published  by  Kaklyugln  and  Norman 
(Reference  3*35)  who,  starting  from  a  Lorentz  model  with  a 
screened  potential,  develop  an  approach  that  accounts  for  the 
spatial  correlation  of  the  colliding  charged  particles.  This 
latter  fact,  not  accountlble  in  the  Boltzman  transport 
equation  approach,  leads  to  a  reduction  of  the  calculated 
conductivity.  This  is  a  result  of  the  localization  of  some 
of  the  free  electrons  that  is  accounted  for  in  this  approach. 
Within  the  context  of  a  model  plasma  (singly-ionized,  non¬ 
degenerate)  they  are  able  to  devise  limiting  scaling  laws 
that  express  the  ratio  as  a  function  of  the  number  of 

particles  (  A/p  )  in  the  Debye  sphere  (  «■,  is  the  ideal-gas 
conductivity  valid  if  A/p  »  /  ).  This  is  clearly  a  valuable 
piece  of  work,  though  as  published  it  assumes  a  somewhat 
idealized  system  as  compared  to  the  real  situation  namely: 
neglect  of  degeneracy,  assumption  of  a  singly  ionized  plasma 
etc . 
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The  results  of  recent  Soviet  experimental  efforts  to 
determine  electrical  conductivity  present  some  ambiguities 
that  need  to  be  resolved.  Batenin  and  Minaev  (Reference  3-31) 
conclude  on  the  basis  of  arc  measurements  In  argon  and 
mercury  that  modified  versions  of  the  Spltzer  theory  (Refer¬ 
ence  3.41)  which  are  supposed  to  be  valid  down  to  A *4-6, 
show  serious  discrepancies  from  the  measured  results  starting 
at  /»iA=4  (  A*  55).  On  the  other  hand  Mltln  (Reference 
3.38)  working  with  argon  and  copper  In  pulsed  arcs  find 
Spltzer  conductivity  results  to  be  adequate  even  for 
7^A=1.6.  Mltln  also  suggests  an  Interesting  semi-emplrlcal 
modification  of  Spltzer  conductivity  formula  based  on  the 
following  changes:  (1)  replace  A  by  (i  +  A*/^,  and  (2) 
define  A~  n^^/b.  In  place  of  the  usual  A=  D/k  for  cases 
where  >  D. 

3.6  Conclusions  of  Theoretical  Program 

The  approximate  values  of  characteristic  plasma  parameters 
have  been  established  by  applying  Ideal  plasma  theory  and  the 
Debye-Huckel  theory  of  weakly  nonldeal  plasma  to  a  nominal 
explosive  driver  condition:  10  kbar  xenon  driver  with  8.5 
km/sec  piston  velocity.  Prom  a  review  of  the  recent  litera¬ 
ture  on  nonldeal  plasmas  together  with  these  preliminary 
calculations  we  have  come  to  the  following  conclusions: 
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1.  Plasma  produced  by  an  explosive  driver  Is  moderately 
to  strongly  nonideal,  with  Interaction  parameter  /  of  order 
unity . 


2.  The  thermodynamics  and  electrical  conductivity  of 
moderately  and  strongly  nonideal  plasmas  are  not  well  under¬ 
stood  at  present. 

3.  The  theories  applicable  to  Ideal  and  weakly  nonideal 
plasmas  can  not  be  extrapolated  to  treat  the  moderately  and 
strongly  nonideal  gases.  Quantum  effects  play  a  dominant  role 
In  such  regimes  and  only  theories  that  take  adequate  account 
of  these  effects  (diffraction  and  exchange  effects  In  close- 
range  scattering  and  electron  degeneracy  effects)  will  give 
meaningful  results.  This  Is  particularly  true  for  strongly 
nonideal  plasmas  where  quantum  effects  are  needed  for  guaran¬ 
teeing  thermodynamic  stability  and  proper  treatment  of  phase- 
transition  phenomena.  Including  plasma  metallization. 

U.  The  Soviets  are  currently  conducting  an  aggressive 
and  wide-ranging  effort  to  resolve  these  problems  in  nonideal 
plasmas.  Their  approach  covers  three  major  areas:  (1)  a 
strong  effort  In  basic  theoretical  formulation  and  numerical 
computation;  (2)  a  variety  of  experiments  to  obtain  valid 
data  in  nonideal  plasma  regimes;  and  (3)  the  formulation  of 
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seml-emplrlcal  techniques  to  provide  useful  predictive  methods 
In  nonideal  plasmas. 

5.  Most  of  the  published  Soviet  work  has  been  based  on 
a  model  plasma  —  a  singly  Ionized,  nondegenerate  system  — 
and  their  results  must  be  extended  to  Include  multiple 
Ionization  and  Incipient  degeneracy  effects  for  explosively 
driven  plasmas. 

3.7  Preliminary  Estimates  of  Explosive  Plasma  Conductivity 

Using  estimated  thermodynamic  properties  reported  In  pre¬ 
vious  sections,  we  have  made  preliminary  estimates  of  the 
electrical  conductivity  of  xenon  at  a  nominal  explosive  driver 
condition  of  8.5  km/sec  flow  velocity  and  10  kbar  pressure. 

For  these  calculations  we  have  used  three  techniques  that  are 
applicable  to  multiply-ionized  plasmas  and  enable  rapid  com¬ 
putation  of  the  conductivity  from  the  simple  analytical 
expressions  given  below: 

A.  Spltzer-Harm  Theory  (Reference  3.24) 

o-  «J2S. <*  ( (mho /cm  )  (3.21) 

where  a  Is  a  tabulated  function  depending  on  the  average 
Ionic  charge  z  (•  n/  ZZ^nf)  and  7  Is  expressed  in  e/ . 

^  4 
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B.  Spltzer-Harm-Mltln  Theory  (Reference  3*38) 

Same  as  above  with  following  modifications:  a) 

replaced  by  ;  and  b)  average  charged  particle 

separation  ''c  exceeds  the  Debye  length  Dc  then  the  plasma 
parameter  Is  redefined;  place  of  the  usual 

definition  A.  - 

C.  High  Temperature  Limit  of  Rogov’s  Theory  (Reference 

3.30) 


a-  =  3£.9r  /3  (z)/  2  ji>7  (3-22) 

Where  2  Dc  .  /9  ,  like  a  In  Spltzer's  expression 

above.  Is  a  tabulated  function  of  the  average  Ionic  charge 

z  . 

Values  of  conductivities  calculated  using  these  three 
methods  are  shown  In  Table  3.6.  Thermodynamic  properties 
evaluated  using  both  the  Ideal-gas  assumption  and  the  Debye- 
Huckel  theory  without  Ionic  shielding  have  been  used  for 
each  conductivity  technique.  The  results  are  seen  to  be 
strongly  dependent  on  the  particular  conductivity  theory 
chosen  and  especially  noteworthy  Is  the  order  of  magnitude 
discrepancy  exhibited  between  the  prediction  based  on  the 
Spltzer-Harm  theory  and  that  of  Rogov's  limiting  expression 
(for  the  Debye-Huckel  thermodynamics  case). 
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TABLE  3.6 

PRELIMIMARY  ESTIMATES  OP  ELECTRICAL  CONDUCTIVITY 
10  KBAR  XENON  DRIVER 


Ideal  Plasma 

Debye-Huckel 
without  Ion 

Approximation 

Shielding 

Spltzer-Harm 

1920  mho/cm 

4886  mho/cm 

Spltzer-Karm-Mltln 

1439 

1237 

Rogov 

529 

435 
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Keeping  In  mind  that  experimental  measurements  of 
conductivity  for  this  particular  plasma  condition  (xenon, 

10  kllobar,  8.5  km/sec)  in  an  imploding  cylinder  device 
yielded  values  In  the  range  of  300-3^(0  mho/cm,  we  can  make 
the  following  observations: 

(1)  The  Spltzer-Harm  theory  Is  totally  inapplicable 
to  the  calculation  of  electrical  conductivity  at  typical 
explosive  driver  plasma  conditions.  This  conclusion  is  not 
surprising  In  view  of  the  fact  that  the  Spltzer-Harm  theory 
assumes  a  classical  model  for  charged  particle  collisions 
and  furthermore  represents  an  asymptotic  solution  valid  for 

A  >'>  \  • 

(2)  Ad  hoc  modifications  of  the  Spltzer-Harm  theory 
such  as  those  proposed  by  Mltin  are  likely  to  be  only  of 
limited  usefulness  at  plasma  conditions  of  Interest  to  us. 

(3)  Rogov’s  theory  shows  considerable  promise  as  a 
rapid  method  for  estimating  conductivities  at  explosive 
driver  conditions.  This  approach  Is  based  on  a  fourth- 
order  Chapman-Enskog  solution  of  the  Boltzmann  equation. 
However,  in  place  of  a  classical  model  for  charged-particle 
collisions  used  In  the  Spltzer-Harm  theory  and  the  various 
modifications  applied  to  it  (Reference  3*^1) >  Rogov  uses  an 
approximate  quantum-mechanical  collision  model  based  on  the 


f] 

U 
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Born  Approximation.  The  striking  improvement  in  agreement 
with  experimental  measurements  exhibited  by  Rogov's  theory 
clearly  establishes  the  Importance  of  accounting  for  o.uantum 
mechanical  colllslonal  effects  in  the  theoretical  description 
of  the  plasma. 

Rogov's  theory  is,  however,  not  rigorous.  Firstly, 
it  treats  the  quantum  colllslonal  effects  in  only  an  approxi¬ 
mate  fashion.  Secondly,  because  it  is  based  on  the  Boltzmann 
transport  equation  approach,  it  falls  to  allow  for  various 
important  effects  associated  with  the  nonideal  behavior  of 
the  plasma.  For  example  it  neglects  the  possibility  of 
spatial  correlation  of  the  charged  particles  which,  as  pointed 
out  by  Kaklyugln  and  Norman  (Reference  3.35)  may  lead  to 
localization  of  some  of  the  free  electrons  in  a  nonideal 
plasma  and  thereby  tend  to  lower  the  conductivity.  In  view 
of  these  limitations,  it  is  clearly  necessary  to  make  a 
careful  study  of  the  limits  of  applicability  of  the  Rogov 
theory.  There  is  undoubtedly  considerable  scope  for  improve¬ 
ment  of  this  theory  to  the  point  where  it  could  be  relied 
upon  as  an  accurate  predictor  of  conductivity  over  most  if 
not  all  of  the  plasma  parametric  regime  of  Interest  in 
explosive  driver  development.  Clearly,  concurrent  experimental 
efforts  toward  obtaining  accurate  conductivity  measurements 
would  greatly  aid  the  achievement  of  this  important  goal. 
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4.  CONCLUSIONS  AND  RECOMMENDATIONS 

4.1  Summary  of  Program  Accomplishments 

1.  Experiments  were  successfully  performed  using  three 
explosive  driver  types  (Imploding  cylindrical,  exploding 
cylindrical,  and  disc)  and  two  gases  (argon  and  xenon). 

2.  Plasma  diagnostic  data  was  successfully  obtained 
using  an  MHD  velocity  gage,  an  eddy  current  plasma  conduc¬ 
tivity  gage,  plasma  resistance  measurements,  an  optical 
pyrometer,  and  shock  wave  instrumentation. 

3.  A  literature  survey  on  nonldeal  plasma  theory  and 
data  was  carried  out,  and  considerable  insight  was  gained 
regarding  the  physics  of  explosively  driven  plasmas. 

4.  Preliminary  calculations  of  plasma  conductivity 
using  several  theoretical  models  were  carried  out  and  com¬ 
pared  with  experimental  data.  The  results  demonstrate  the 
importance  of  nonldeal  plasma  phenomena,  and  provide 
encouraging  progress  in  removing  the  discrepancy  between 
experiment  and  theory . 

4.2  Recommendations 

1.  Further  experimental  work  needs  to  be  performed  to 
provide  an  adequate  data  base  on  electrical  conductivity, 
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temperature,  and  other  accessible  parameters  of  explosively 
driven  plasmas. 

2.  Experimental  and  theoretical  work  should  be  per¬ 
formed  on  ablation  and  boundary  layer  growth  In  explosively 
driven  plasmas,  with  particular  attention  paid  to  electrical 
conductivity  and  velocity  gradients  In  the  plasma. 

3.  To  provide  accurate  experimental  data  further 
diagnostic  development  work  must  be  performed.  In  particular 
three-dimensional  MHD  calculations  should  be  performed  to 
Investigate  the  conductivity  averaging  process  In  the  eddy 
current  conductivity  gage  and  plasma  resistance  measurement. 
Additional  effort  to  reduce  nonideal  effects  In  the  explosive 
driver  would  also  Improve  accuracy.  Some  development  of  the 
optical  pyrometer  Is  necessary  to  obtain  absolutely  cali¬ 
brated  data. 

4 .  Further  work  Is  required  to  develop  the  theory  of 
nonideal  plasmas  to  the  point  where  useful  predictions  of 
experimental  environments  can  be  made. 


-118- 


FR-119 


REFERENCES 


1.1  Gill,  S.  P.  et  al,  "Explosive  MHD  Research",  Artec 
Associates  Inc.  Annual  Report  119AR,  Navy  Contract 
N0001il-75-C-0822,  April  1976. 

1.2  Baum,  D.  W.  et  al,  "Development  of  High  Energy  Density 
Simulator",  Artec  Associates  Inc.  Final  Report  120, 
Defense  Nuclear  Agency  Contract  DNA001-75-C-0271 , 
December  1976. 

1.3  Baum,  D.  W.  et  al,  "Explosive  MHD  Power  Source",  Artec 
Associates  Inc.  Final  Report  113,  Navy  Contract  N00178- 
75-C-0120,  June  1977. 


3.1  Zeldovlch,  Ya.  B  and  Ralzer,  Yu.  P.,  "Physics  of 

Shock  Waves  and  High  Temperature  Hydrodynamic  Phenomena", 
Academic  Press,  New  York,  1966. 

3.2  Berlin,  T.  H.  and  Montroll,  E.  W. ,  "On  the  Free  Energy 
of  a  Mixture  of  Ions:  An  Extension  of  Kramer's  Theory", 
J.  Chem.  Phys.,  1952,  p.  75. 

3.3  Grlem,  H.  R. ,  "High  Density  Corrections  In  Plasma 
Spectroscopy",  Phys.  Rev.,  128,  1962,  p.  997. 

3.^  McChesney,  M.,  "Equilibrium  Shockwave  Calculations  In 
Inert  Gas,  Multiply-Ionized  Debye-Huckel  Plasmas", 

Can.  J.  Phys.,  196*»,  p.  2U73. 

3.5  Vedenov,  A.  A.,  "Thermodynamics  of  a  Plasma",  Reviews 
of  Plasma  Physics,  Vol.  1,  Consultants  Bureau,  New 
York,  1965,  p.  312. 

3.6  Stewart,  J.  C.  and  Pyatt,  K.  D.,  "Lowering  of  the 
Ionization  Potential  In  Plasmas",  Astrophys.  J.  1Jl4, 

1966,  p.  1203. 

3.7  Van  Thiel,  M.  and  Alder,  B.  J.,  "Shock  Compression  of 
Argon",  J.  Chem.  Phys.,  1966,  p.  IO56. 

3.8  Brush,  S.  G.  et  al ,  "Monte  Carlo  Study  of  a  One- 
Component  Plasma",  J.  Chem.  Phys.,  1966,  p.  2102. 

3.9  Norman,  G.  E.  and  Starostln,  A.  N.,  "Insufficiency 

of  the  Classical  Description  of  a  Nondegenerate  Dense 
Plasma",  Tepl.  Vys.  Temp.,  6,  1968,  p.  410. 


-119- 


PR-119 


3.10  Ross,  M. ,  "Shock  Compression  of  Argon  and  Xenon  TV: 
Conversion  of  Xenon  to  a  Metal-like  State",  Phys.  Rev., 
171,  1968,  p.  777. 

3.11  Plchakchi,  L.  D. ,  "The  Thermodynamics  of  a  Dense 
Plasma",  Sov.  Phys.  JETP,  2^,  1968,  p.  8^5. 

3.12  Rakhimov,  A.  T.  and  Ullnlch,  P.  R.,  "The  Possibility 

of  Exlstense  of  a  Strongly  Nonldeal  Plasma",  Sov.  Phys. 
JETP,  30,  1968,  p.  772. 

3.13  Vorob’ev,  V.  W.,  et  al,  "Statistical  Theory  of  a  Dense 
Plasma",  Sov.  Phys.  JETP,  30,  1970,  p.  459. 

3.14  Norman,  G.  E.  and  Starostin,  A.  N.,  "Thermodynamics  of 
a  Strongly  Nonldeal  Plasma",  Tepl.  Vys.  Temp., 

1970,  p.  413. 

3.15  Ermokhln,  N.  V.  et  al ,  "An  Experimental  Investigation 
of  Coulomb  Interaction  In  a  Dense  Plasma",  Tepl.  Vys. 
Temp.,  9,  1971,  p.  665. 

3.16  Lomakin,  B.  N.  and  Portov,  V.  E. ,  "Pulsed  X-ray  Studies 
of  Shockwaves  in  Dense  Cesium  Vapor",  Tepl.  Vvs.  Temp., 
9,  1971,  p.  1291. 

3.17  Portov,  V.  E.  et  al,  "Thermodynamic  Properties  of  a 
Cesium  Plasma",  Tepl.  Vys.  Temp.,  9,  1971,  p.  869. 

3.18  Lomakin,  B.  N.  and  Portov,  V.  E.,  Zh.  Eksp.  Fix.,  63 , 
No.  7.  1972. 


3.19  Zelener,  B.  V.  et  al,  "Statistical  Theory  of  a  Nonldeal 
Plasma",  Tepl.  Vys.  Temp.,  1972,  p.  II60. 

3.20  Rogers,  F.  J.  and  Dewitt,  H.  E. ,  "Statistical  Mechanics 
of  Reacting  Coulomb  Gases",  Phys.  Rev.  A,  1973, 

p.  1061. 


3.21  Zelener,  B.  V.  et  al,  "Thermodynamics  and  Correlation 
Function  of  Weakly  Nonldeal  Non-Debye  Coulomb  Systems", 
Tepl.  Vys.  Temp.,  U,  1973,  p.  922. 

3.22  Zelener,  B.  V.  et  al,  "Equations  of  State  and  Ionization 
Nonequilibrium  of  a  Nonldeal  Plasma",  Tepl.  Vys.  Temp., 
13,  1975,  p.  913. 


3.23 


Zelener,  B.  V.  et  al,  "Pseudopotential  Model  of  a 
Nonldeal  Plasma",  Temp.  Vys.  Temp.,  13,  1975,  p.  712. 


-120- 


FR-119 


3.2il  Spitzer,  L.,  "Physics  of  Fully  Ionized  Gases",  Wiley 
Interscience,  2nd  Edition,  1962. 

3.25  Ahtye,  W.  P.,  "A  Critical  Evaluation  of  Methods  for 
Calculating  Transport  Coefficients  of  Partially  and 
Fully  Ionized  Gases",  NASA  TN-D  2611,  January  1965* 

3.26  Mason,  E.  A.  et  al ,  "Transport  Coefficient  of  Ionized 
Gases",  Phys.  Fluid,  1967,  p.  1327. 

3.27  Robinson,  J.  V.,  "Measurement  of  Plasma  Energy  Density 
and  Conductivity  from  3  to  120  kbar" ,  J.  Appl.  Phys. 

38,  1967,  P-  210. 

3.29  Williams,  R.  H.  and  Dewitt,  H.  W.,  "Quantum  Mechanical 

Plasma  Transport  Theory",  Phys.  Fluid,  1^,  1969,  p.  2326. 

3.29  Bakeev,  A.  A.  and  Rovlnsky,  R.  E.,  "Electrical  Pro¬ 
perties  of  High-Pressure  Pulsed  Discharge  Plasma  In 
Inert  Gases",  Tepl .  Vys.  Temp.,  8,  1970,  p.  1121. 

3.30  Rogoff,  V.  S.,  "Calculation  of  Plasma  Conductivity", 

Tepl.  Vys.  Temp.,  8,  1970,  p.  689. 

3.31  Batenin,  V.  M.  and  Minaev,  D.  V.,  Tepl.  Vys.  Temp., 

9,  1971,  p.  676. 

3.32  Povlnskli,  R.  E.,  "Electrical  Conductivity  of  a  Non- 
Debye  Plasma",  Tepl.  Vys.  Temp.,  1972,  p.  1. 

3-33  Barolskll,  S.  G.  et  al,  "Measurement  of  the  Electrical 
Conductivity  of  a  Dense  Strongly  Nonideal  Cesium 
Plasma",  Sov.  Phys.  JETP,  35,  1972,  p.  9^. 

3.3^  Devoto,  R.  S.,  "Transport  Coefficients  of  Ionized 
Argon",  Phys.  Fluid,  16,  1973,  p.  616. 

3.35  Kaklyugin,  A.  S.  and  Norman,  G.  E.,  "Electrical  Con¬ 
ductivity  of  a  Non-Debye  Plasma",  Tepl.  Vys.,  Temp., 

11,  1973,  p.  238. 

3.36  Litvinov,  I.  I.,  "Electrical  Conductivity  of  a  Non- 
Debye  Plasma",  Tepl.  Vys.  Temp.,  3^,  197'^,  p.  1100. 

3.37  Andreev,  S.  and  Gavrilova,  T.  V.,  "Measurement  of 
Electrical  Conductivity  of  Air  Plasma  at  Pressures 
Above  100  atm.".  Tempi.  Vys.  Temp.,  1975,  p.  176. 


-121- 


3 -35  Mltln,  R.  B.,  "Experimental  Determination  of  the 
Electrical  Conductivity  of  Dense  Plasmas  with  Low 
*  Coulomb  Logarithm",  Tepl .  Vys.  Temp.,  5^,  1975,  p.  706. 


3.39  Ecker,  G.  and  Kroll,  W.,  "Lowering  of  the  Ionization 
Potential  for  a  Plasma  In  Thermodynamic  Equilibrium", 
Phys.  Fluid,  6,  1962,  p.  62. 

3.^0  Conwell,  E.  M.  and  Welsskopf,  V.  F.,  Phys.  Rev.,  77 
1950,  p.  388.  ~ 

3.^1  Kruger,  C.  H.  et  al,  "Transport  Properties  of  MHD 
Generator  Plasmas",  AIAA  J.,  6,  1968,  p.  1712. 


i 

I 


[ 

I 


